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ABSTRACT 

This  thesis  is  a  study  of  quasi-static,  axially  symmetric 
plastic  deformation  of  a  rigid-plastic,  non-hardening  material  which 
obeys  the  Tresca  yield  criterion.  Investigations  of  the  various 
plastic  regimes  possible  as  the  stress  point  traverses  the  Tresca 
yield  locus  are  undertaken  and  it  is  shown  that,  in  all  non-trivial 
cases,  the  equations  governing  the  associated  stress  and  velocity 
fields  of  each  plastic  regime  are  hyperbolic.  Particular  attention 
is  given  to  the  regime  representative  of  the  Haar-Karman  hypothesis. 
In  view  of  the  success  of  this  hypothesis  in  the  solution  of  many 
problems  in  axially  symmetric  indentation,  it  has  been  used  in  the 
thesis  with  regard  to  an  indentation  problem  with  a  cone. 

A  rigid-plastic,  non-hardening  material,  which  is  semi¬ 
infinite  in  extent,  is  indented  by  rigid,  smooth,  right  circular 
conical  punch  initially  inserted  into  a  prepared  cavity  at  the 
surface  of  the  material.  By  using  the  plastic  stress  field  of  the 
"incomplete"  solution,  a  value  of  4.6424nkR2  is  obtained  for  the 

I 

yield  point  load,  where  k  is  the  maximum  shearing  stress  for  the 
material  and  R  is  the  surface  radius  of  the  prepared  cavity. 


(il) 


ACKNOWLEDGEMENTS 

The  author  wishes  to  extend  his  appreciation  to  Dr. 
J.  B.  Haddow  for  his  guidance  and  supervision  of  this  thesis. 

Additional  thanks  are  extended  also  to  those  members 
of  the  University  of  Alberta  Computing  Centre  whose  assistance 
proved  quite  helpful. 


TABLE  OF  CONTENTS 


Page 

ABSTRACT  oo.oooo.o0„0..  (0 

ACKNOWLEDGEMENTS  ..................  (ii) 

CHAPTER 

I„  INTRODUCTION 

1.1  Historical  Survey  .............  1 

1.2  Scope  of  Thesis  ...  ..........  .  3 

II.  GENERAL  FIELD  EQUATIONS  FOR  PLASTIC  FLOW  UNDER  QUASI- 
STATIC  AXIALLY  SYMMETRIC  CONDITIONS 

2.1  Conditions  of  Axial  Symmetry  .......  5 

2.2  Physical  Conditions  of  Yield  and  Plastic  Flow  8 

III.  ANALYSIS  OF  THE  FIELD  EQUATIONS 

3.1  Group  1 1  Plastic  Regimes  B  and  E .  14 

(a)  Stress  Fields  ..............  14 

(b)  Velocity  Fields  .............  15 

3.2  GroupXIs  Plastic  Regimes  AB  and  EF  .  .  .  .  15 

(a)  Velocity  Field  for  AB  .........  .  16 

(b)  Stress  Field  for  AB  ..........  .  21 

3.3  Group  Ills  Plastic  Regimes  A  and  F  .  .  .  .  24 

(a)  Stress  Field  for  F  ..........  .  24 

(b)  Velocity  Field  for  F  .........  .  27 

3.4  Group  IV s  Plastic  Regime  AF  ......  .  3 0 

(a)  Velocity  Field  for  AF..........  J>0 

% 

(b)  Stress  Field  for  AF  .  .  .  .  . .  3* 


CHAPTER 

IV.  INCIPIENT  PLASTIC  FLOW  IN  A  SEMI- INFINITE  REGION  OF  RIGID- 
PLASTIC  NON -HARDEN ING  ISOTROPIC  MATERIAL  DUE  TO  A  LOAD 
APPLIED  BY  A  SMOOTH  CONICAL  XND ENTER  FITTED  INTO  A  CONICAL 
CAVITY o 

4.1  Preliminary  „  .  0  „  .  „  .  „  0  „  .  „  „  0  .  33 

4.2  Formulation  of  the  Problem  .  .  .  0  .  „  .  .  35 

4.3  Plastic  Stress  Field  ABC  .  .  .  .  .  .  .  .  .  37 

bob  Plastic  Stress  Field  ACD  .........  44 

4.5  Plastic  Stress  Field  ADE  .  .  .  .  .  .  .  .  .  48 

4o6  Calculation  of  Indentation  Stress  from  Slip-Line 

Field  of  the  Incomplete  Solution  .....  50 

4*7  Discussion  and  Future  Requirements  ....  52 

APPENDIX 

I.  COMPUTER  PROGRAMS  AND  OUTPUT  DATA  FOR  PLASTIC  STRESS  FIELDS 
A.  Plastic  Stress  Field  AC^D  .........  5b 

1.  Procedure  .................  5b 

2o  Fortran  Source  Program,  Input  and  Output  Data  55 

Bo  Plastic  Stress  Field  ADE  .  .  .  .  .  .  .  .  .  .  56 

f25 

II.  THE  EVALUATION  OF  /  cr^rdr  BY  NUMERICAL  METHODS  80 


BIBLIOGRAPHY 


85 


, 


LIST  OF  TABLES 


Table 


Page 


I.  Yield  Conditions  and  Flow  Rule  for 

Individual  Plastic  Regimes  . 11 


II.  Values  of 


a 


R 

ao 

Pa, 


and 


R 


*0 


po 


*  A,  I 

III.  Values  of  — f-  and  cpAi 


,46 


LIST  OF  FIGURES 


FIGURE 


1„  Cylindrical  Polar  Coordinate  System  and  Stress 
Components  for  Axially  Symmetric  Conditions  „  , 


Page 


2,  The  Tresca  Yield  Condition 


o  o  o 


3«  Slip~Line  Field  in  a  Meridian  Plane  for  Indentation 
by  a  Smooth  Conical  Punch  Fitted  in  a  Prepared 


Cavity 


00000000000000000000900 


36 


l+o  Reference  Diagram  Illustrating  Procedure  Used  for 

Determination  of  Slip°Line  Field  of  Figure  3  .  .  .  .  38 


5.  Straight-Line  Free  Boundary 


1+0 


dR 


5a,  Diagram  for  Determining  -7-^ 

(XV/ 


dR 


5b,  Diagram  for  Determining 


dy 


00000000000 


00000000 


1+0 

1+0 


60  Intersection  of  Characteristics  Through  Neighbouring 


Points  to  Show  Approximation  R^  for  True  Inter- 


Section  R„ 


0000000000000000000 


1+7 


7,  Reference  Diagram  for  Calculation  of  Indentation 


Stress 


50 


1 


CHAPTER  X 
INTRODUCTION 


1.1  Historical  Survey. 

The  work  of  Tresca  on  the  flow  of  metals  during  punching  and 
extrusion  is  considered  the  beginning  of  the  mathematical  theory  of 
plasticity.  Tresca  [1864]  advanced  the  theory  that  a  metal  yields 
plastically  when  the  maximum  shearing  stress  attains  a  critical  value. 

This  hypothesis  is  now  called  the  Tresca  yield  criterion.  Saint-Venant 
[1870]  used  this  criterion  to  interpret  the  results  of  experiments  on 
the  torsion  and  bending  of  cylinders.  Saint-Venant  proposed  that  during 
two  dimensional  plastic  flow  in  isotropic  materials  a  co-axial  relationship 
exists  between  the  stress  tensor  and  the  strain-rate  tensor.  Levy  [I87O] 
extended  this  idea  of  Saint-Venant 5 s  to  three  dimensions  and  postulated 
proportionality  of  the  stress  deviator  components  and  the  corresponding 
plastic  strain-rate  components.  Von  Mises  [I9I3]  suggested  the  same 
proportionality  relationship  (not  knowing  of  Levy's  earlier  work)  and 
postulated  a  new  criterion  for  the  yielding  of  metals.  The  Mises  yield 
criterion  implies  that  yielding  of  a  material  begins  when  the  octahedral 
shearing  stress  reaches  a  critical  value.  Von  Mises  extended  his  work  in 
1928  to  perfectly  plastic  solids  having  a  general  yield  function.  He 
derived  the  plastic  stress-strain  relations  corresponding  to  this  yield 
function  and  from  this  deduced  the  concept  of  the  plastic  potential  and  its 
associated  flow  rule.  Through  the  independent  works  of  Melan  [1938]  and 
Prager  [  19^4-9 3 »  the  general  plastic  stress-strain  relations  for  solids  with 
a  regular  yield  function  were  obtained.  Meanwhile*  Reuss  [1933]  ^acl 
investigated  the  flow  rule  associated  with  the  Tresca  singular  yield 
criterion.  Koiter  [1953a]  adapted  the  theory  of  plastic  potential  to 
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materials  with  singular  yield  fusieticras  and  this  resulted  in  the  Koiter- 
Prager  generalization  of  the  Mises  theory  of  plastic  potential. 

Application  of  the  mathematical  techniques  associated  with  the 
plastic  potential  has  resulted  in  the  solution  of  many  problems  in  plane 
strain  plasticity.  However,  very  few^non-trivial. axially  symmetric  problems 
in  plasticity  have  been  solved.  Investigators  such  as  Hill  [1948], 

Syraonds  [194-9]  and  Parsons  [1956]  have  shown  that  when  the  Mises  yield 
criterion  and  its  associated  flow  rule  are  used  in  axially  symmetric 
problems,  the  plastic  stress  and  velocity  equations  are  not  hyperbolic. 
Serious  mathematical  difficulties  arise  in  finding  solutions  to  the 
problems  under  such  conditions.  Koiter  [1955^3  the  Tresca  yield 

criterion  and  associated  flow  rule  to  obtain  closed  solutions  for  many 
problems  of  partially  plastic,  thick-walled  tubes  under  axial  end-conditions. 
Shield  [19553  solved  the  axially  symmetric  problem  of  incipient  plastic 
flow  of  a  semi-infinite,  non-hardening,  rigid-plastic  material  indented  by 

a  rigid,  smooth^cy lindrical  indenter.  In  his  analysis,  Shield  used  the 

.  1 

Tresca  yield  criterion  and  its  associated  flow  rule  together  with  the  Haar- 

von  Karman  hypothesis.  This  hypothesis  had  been  under  criticism  by  many 

authorities  of  plasticity;  e.g^Hill  [1950a],  since  its  inception  in  1909. 

Haar  and  von  Karman  [19093  postulated  that  in  some  statically  determinate 

problems  of  axial  symmetry,  the  circumferential  stress  is  equal  to  one  of 

the  principal  stresses  in  the  axial  plane.  Application  of  this  hypothesis 

without  justification  was  the  main  criticism.  Ishlinskii  [ I9I4.I4-]  attempted 

the  problem  of  indentation  of  the  plane  surface,  of  a  semi-infinite  material 

by  a  circular,  flat-ended,  smooth,  rigid  punch.  Berezancev  [1955]  attempted 

the  problem  in  soil  mechanics  of  normal  penetration  of  cohesive  soils  by  a 

rigid,  smooth,  right  circular  cone.  Both  Ishlinskii  and  Berezancev  assumed 
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the  Haar-von  Karman  hypothesis.  The  criticism  directed  to  Ishlinskii’s 
work,  other  than  the  ac  hoc  use  of  the  Haar-von  Karman  hypothesis,  was  the 
inherent  inaccuracy  of  the  graphical  method  used  in  obtaining  the  plastic 
stress  field  and  in  not  attempting  to  find  the  associated  plastic  velocity 
field.  Berezancev"s  work  is  only  approximate  due  to  his  assumed  conditions 
on  the  boundary.  Also,  there  was  no  attempt  made  to  find  an  associated 
velocity  field.  It  was  not  until  Shield  [1955]  had  derived  the  exact 
solution  of  Ishlinskii y s  problem  that  these  works  were  considered  justified 
in  utilizing  the  Haar-von  Karman  hypothesis. 

1.2  Scope  of  Thesis. 

In  this  thesis,  the  Koiter-Prager  generalization  of  the  von  Mises 
theory  of  plastic  potential  is  used  to  derive  the  quasi-static  stress  and 
velocity  equations  for  a  material  deforming  plastically  under  conditions 
of  axial  symmetry.  The  material  obeys  the  Tresca  yield  criterion  and  is 
assumed  to  be  rigid-plastic,  non-hardening  and  isotropic  in  nature.  The 
various  plastic  regimes  possible  as  the  stress  point  traverses  the  Tresca 
yield  locus  are  conveniently  represented  as  members  of  four  distinct  groups. 
In  each  group,  the  field  equations  for  the  associated  plastic  regimes  are 
shown  to  be  non-elliptic . 

The  hyperbolic  stress  equations  of  the  plastic  regime  represented 
by  the  Haar-von  Karman  hypothesis  are  used  in  an  indentation  problem  with 
a  cone.  A  material,  semi- inf inite  in  extent,  rigid-plastic  and  non¬ 
hardening,  has  on  its  stress  free  plane  surface  a  right  conical  cavity 
whose  axis  is  normal  to  the  plane,  surface.  A  rigid,  smooth,  right  circular 
conical  indenter  is  inserted  into  this  cavity  occupying  it  fully.  The 
indenter  is  then  normally  loaded  until  plastic  flow  of  the  material  occurs. 


-  k  - 

The  plastic  stress  field  and  the  yield  point  load  for  incipient  plastic 
flow  of  the  material  is  determined  by  methods  of  numerical  analysis.  No 
attempt  was  made  to  derive  an  associated  kinematically  admissible  velocity 
field  or  to  extend  the  stress  field  into  the  remaining  rigid  region. 
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CHAPTER  II 

GENERAL  FIELD  EQUATIONS  FOR  PLASTIC  FLOW  UNDER 
QUASI-STATIC  AXIALLY  SYMMETRIC  CONDITIONS. 


2.1  Conditions  of  Axial  Symmetry. 

Let  0  be  the  origin  of  a  right-handed  3"dimensional  system  of 
cylindrical  polar  coordinates  r,  Q ,  z  (Fig.  l).  With  respect  to  this 


Cylindrical  Polar  Coordinate 
System  and  Stress  Components 
for  Axially  Symmetric 
Conditions  0 


system,  the  physical  components  of 
the  stress  tensor  are  denoted  as 


(cr  ,  or  ,  a  ,  cr.  ,  cr  ,  cr  _),  the 
v  r*  Q  z  Qz  rz  r0'’ 

physical  components  of  the  strain- 
rate  tensor  as 


(e,€^,e,€^,e  ,€^),  and 

v  r  9  z  Qz  rz  r G 

the  physical  components  of  the 
velocity  as  (u,  v,  w) . 


Under  conditions  of  axial  symmetry,  the  choice  of  reference  plane 

is  arbitrarily  any  meridian  plane.  The  shear  components  <j  and  <j  , 

Qz  r  9 

the  velocity  component  v  ,  and  the  shear  strain-rate  components  e  and 

QZ 

e  must  be  zero.  The  remaining  components  of  stress,  velocity  and  strain- 
r0 

rate  are  then  expressible  as  functions  of  r,  z  and  T  where  T  is 
physical  time.  However,  in  any  quasi-static  problem,  a  time-scale  is 
necessary  only  to  order  events  in  contrast  to  dynamical  problems  where 
inertial  effects  are  considered  and  the  physical  time  must  be  used.  Thus 
in  quasi-static  problems,,  any  suitable  monotonically  increasing  parameter 
t  correlated  with  progressive  deformation  may  be  used  as  a  time-scale. 
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If  body  forces  are  neglected  and  quasi-static  conditions  are 
assumed,  the  equations  of  equilibrium  satisfied  by  the  stress  components  are 


cr 

r 


r 


0 


(2.1.1) 


and 


(2. 1.2) 


The  physical  components  of  strain-rate  are  expressible  as 


e  = 


ehj 

dr 


e  = 

e 


u 


dw 

dz 


rz 


(2.1 


Along  the  axis,  r  =  0  ,  conditions  are  imposed  upon  the  components 

of  stress,  velocity,  strain-rate  and  their  respective  derivatives.  Derivation 

of  these  conditions  are  based  on  the  assumption  that  plastic  flow  occurs 

without  fracture.  Also? conditions  must  be  such  as  to  insure  the  existence 

of  derivatives  and  certain  limits  of  stress,  velocity  and  strain-rate  upon 

approaching  the  z-axis.  From  the  equations  of  equilibrium  (2.1.1)  and  (2.1.2), 

it  cannot  be  concluded  that  cr  -  or  =  cr  =0  when  r  =  0  in  order  to 

r  Q  rz 

avoid  infinite  values  of  the  stresses.  Although  stresses  may  be  bounded, 

there  is  no  reason  why  their  derivatives  should  not  become  unbounded  on 

approaching  the  axis.  o*r  ,  at  any  point  P,  is  defined  as  the  normal  stress 

on  a  plane  through  P  perpendicular  to  a  radius  through  P  *  and  cr^  is 

defined  as  the  normal  stress,  at  P,  on  a  plane  through  P  containing  the 

z-axis.  Hence5as  argued  by  Parsons  [1956],  these  definitions  must  hold  for 

all  points  and,  in  particular,  for  those  on  the  z-axis.  When  P  lies  on 

this  axis,  any  plane  through  P  containing  the  axis  is  also  perpendicular 

to  a  radius  drawn  through  P.  Hence  <t  and  cr  at  points  on  the  axis 

r  u 
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are  the  same  stress  and  -  <T^  =  0.  Also  for  any  point  P  on  or  off  the 

z-axis,  or  is  defined  as  the  tangential  stress  at  P  on  a  plane  through 

P  perpendicular  to  the  z-axis.  If  P  lies  on  this  axis,  it  follows  that 

azr  =  0  ,  for  if  a  /  0  this  tangential  stress  on  the  plane  would  have 

a  definite  direction  at  P,  contradicting  the  hypothesis  of  axial  symmetry. 

Hence  conditions  of  axial  symmetry  and  proportionality  of  the  components 
/ 

of  the  stress  and  plastic  strain-rate  tensor  require  that 


<7 

zr 


0 


€  = 
r 


e  =  0 
zr 


} 


on  r  =  0  . 


(2.1.4) 


The  condition  that  particles  near  the  axis,  r  =  0  ,  do  not 
separate  requires  that  u  =  0  when  r  =  0. 


From  (2.1.1), 

/  da 


lim 
r  — ►  0 


it  follows  that 


rz 


a  -cr. 

r  e 


0  . 


But  since  cr  =  0  on  r  =  0 
rz 


lim 
r  -*■  0 


^°rz 

=0- Hence 


lim 
r  0 


0 


(2.1.5) 


From  (2.1.2),  it  also  follows  that 


lim 
r  —  0 


lim 

r  -=►  0 


lim 
r  --*•  0 


a 

rz 

r 


Hence 


/ 


3  Oi  OC»  r»« 


lim 
r  0 


dcr 

z 

+ 
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0 


(2.1.6) 


rz 


From  (2o  lo4)  s,  €rz  =  -  +  §7  )  =  0  on  r  =  °' 


Since  u  =  0  on  r  =  0  ,  =  0  on  r  =  0 


6w 

Hence  lim  =  0 

r  — ►  Q  r 


(2.1.7) 


Condition  (2. 1.7)  insures  continuity  of  €^z  at  the  axis  of 
symmetry.  (2. 1 . 15) » (2. 1 . 6) , (2. 1 .7)  are  the  conditions  on  the  limits  of  the 
respective  quantities  which  must  be  satisfied  under  axial  symmetry. 


2.2  Physical  Conditions  of  Yield  and  Plastic  Flow. 

Consideration  is  now  confined  to  a  rigid-plastic,  non-hardening, 
homogeneous  and  isotropic  material  which  obeys  the  Tresca  yield  criterion. 
Rigid-plastic  implies  incompressibility  and  the  vanishing  of  all  elastic 
strains.  Non-hardening  expresses  the  condition  that  non-vanishing  plastic 
strain-rates  may  occur  when  the  material  is  under  a  constant  state  of  stress 
at  the  yield  limit. 


The  Tresca  yield  criterion  of  constant  maximum  shearing  stress 
for  isotropic  materials  is  conveniently  described  in  a  3-dimensional 
principal  stress  space.  In  this  space,  the  principal  stress  components 
(T,  ,  cr0,  <r,  are  chosen  as  rectangular  Cartesian  coordinates  and  any  state 
of  stress  is  represented  by  a  point  o\  (i  =  1,2,3)*  The  Tresca  yield 
criterion  is  represented  by  the  surface  of  a  regular  hexagonal  prism  with 
its  axis  equally  inclined  to  the  positive  CTp  <J^9  or^  axes  and  passing 
through  the  origin.  No  change  in  the  plastic  strains  of  an  element  of  a 
material  occurs  if  the  stress  point  lies  within  the  prism.  Such  states  are 
referred  to  as  "safe".  Increments  of  plastic  strain  can  occur  only  if  the 
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stress  point  is  on  the  surface;  hence,  the  surface  is  called  the  yield 

* 

surface  .  States  of  stress  which  are  either  within  or  on  the  yield  surface 

•  *  •  < 

are  called  "allowable",,  No  state  of  stress  can  lie  outside  the  yield  surface. 
The  intersection  of  this  prism  by  a  general  plane  <j  =  constant  is  an 


Figure  2 

The  Tresca  Yield  Criterion 


irregular  hexagon.  Figure  2  is  the  orthogonal  projection  onto  the  <j ^  ~  cr ^ 
plane  of  such  an  intersection  -  a  distance  cr  from  the  origin  0.  Points 
on  this  irregular  hexagon  ABCDEF  represent  all  possible  states  of  stress  with 


maximum  shearing  stress  k0  Homogeneity  of  the  material  requires  that  any 
element  of  the  material  will  become  plastic  when  the  maximum  shearing  stress 
on  that  element  has  attained  the  value  of  this  material  constant  k. 


According  to  the  concept  of  plastic  potential  (Hill  [1950c]),  the 
principal  plastic  strain~rate  free  vector  2Ge(€^),  in  principal  stress 

space,  associated  with  the  principal  stress  bound  vector  cr(o\)s  has  the 

-  ( 

same  direction  as  the  outward  normal  to  the  regular  yield  surface,  f(o\)  =  0, 
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at  the  point  or  (i  =  1,2,3)  which  represents  the  plastic  state  of  stress. 
The  factor  2G  is  used  only  to  give  dimensions  of  stress  to  the  plastic 
strain-rate  vector  c(€  )  in  the  principal  stress  space.  The  plastic 
strain-rate  components  are  determined  by 


where 


i 


(2,2,1) 


A  =  0  if  f  <  0  and  also  if  f  =  0  and  f  =  v—  <t  <  o  ; 

0(J.  i 
i 

A  >  0  if  f  =  0  and  f  =  0;  A  being  an  indeterminate  scalar 
parameter.  An  upper  dot  associated  with  a  quantity  denotes  differentiation 
with  respect  to  time  or  any  other  monotonically  increasing  parameter  correlated 
with  progressive  deformation. 


Singular  yield  surfaces  with  edges  or  corners  are  represented  by 

a  finite  or  infinite  number  of  yield  surfaces  f  (<x„),  a  =  1,2,,, „  .  States 

a  i 

of  stress  at  the  yield  surface  are  described  by  a  value  zero  of  one  or  more 
of  the  yield  functions^  all  other  yield  functions  being  negative.  By  the 

Koiter-Prager  generalization  of  the  plastic  potential  (Koiter  [i960]),  the 
direction  of  the  plastic  strain-rate  vector  2Ge(e^)  is  unique  and  coin¬ 
cident  with  the  outward  normal  to  the  yield  surface  except  at  singular  points 
along  the  edges.  At  a  singular  point ,  the  direction  of  2Gc(e^)  must  lie 
between,  and  in  the  plane  defined  by,  the  unique  normals  drawn  outwards  to 
the  two  faces  of  the  prism  intersection  at  the  particular  singular  point 
considered.  The  plastic  strain-rate  components  are  determined  by  the  formulae 

df 

e,  =  A “  ,  (i  =  1,2,3)  ,  (2.2.2) 

1  a 

1 


where 


11 


A  =  0  if 

a 

A  >  o  if 
a 

parameters „ 


df 

f  <0  and  also  if  f  =  0S  f  ~  er.  <  0  : 

a  a  a  acr^  i 


f  =0  and  f  =  0;  A  being  indeterminate  scalar 
ct  a  cl 


In  this  thesis,  discussion  of  the  types  of  plastic  flow  possible 
as  the  stress  point  cr^  traverses  the  Tresca  yield  locus  is  confined  to 
EFAB  (Fig.  2)  along  which  >  ov,  .  The  locus  EFAB  is  divided  into  in¬ 
distinct  groups  of  plastic  regimes;  namelys  I,  B  and  E;  II,  AB  and  EF; 
III,  A  and  F;  and  IV,  AF.  Groups  I  and  III  represent  singular  edge 
plastic  regimes  while  groups  II  and  IV  represent  regular  face  plastic 
regimes.  The  yield  functions  corresponding  to  the  regular  plastic  regimes 
are 


"  2k  , 


fg£  -  cr^  ~  +  2k  • 


(2c2o3) 


The  plastic  strain  rates  for  the  k  groups  are  obtained  by  applying  (2„202) 
to  the  yield  functions  (2. 2 <>3)°  They  are  represented  in  Table  I0 


Table  I 

Yield  Conditions  and  Flow  Rules  for 

Group  Plastic  Yield  Condition 

Regime 


Individual  Plastic  Regimes 

Plastic  Strain  Rates 
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°A^ 
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2k,  a  =  ct 
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-x5 

X2 

/ 
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The  circumferential  direction  is  a  principal  direction,  since 

o\-  =  0.  For  purpose  of  discussion* let  ct  ~  a  and  =  ?  .  If  either 

0  z  r  ?  3  0  3  8 

al  ^  °”'2  °r  €1  ^  €2*  then  the  principal  directions  of  stress  and  plastic 
strain-rate  within  the  reference  plane  are  unique  and  coincident  since  the 
material  is  assumed  isotropic.  By  assumption,  cr^  >  <r0  and  by  isotropy, 

€1  “  €2°  ■^ie  positive,  first  and  second,  principal  directions,  cr^  and  cr^ 
are  then  defined  to  make  an  angle  tp,  0  <  Cp  <  jt  with  the  positive  r  and 
z  directions  respectively.  These  principal  stresses  in  the  r  -  z  refer¬ 
ence  plane  are 
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(cr 

+ 
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{f 
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rzJ 
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= 

1 
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)2 

cr2  F 

• 
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'•4 
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rzJ 

(2.2.5) 


The  principal  piastre  strain-rate  components  are 


€.=-(€  +  €  )  +  {r  (e 

1  2  r  z'  l4  v  r 

€0  -  “  ( €  +  €)  =  {r(e 

2  2  r  z '  1 4  '  r 


l 


(2.2.6) 


The  stress  components  in  the  r  and  z  directions  are 


°r  =  2  ^1  +  °2^  +  2  (°1  “  °2^  CO®  ^  * 

s  2  ^1  +  ^  °  2  (°1  “  *2^  C°S  ^  *  (2.2.7) 

arz  *  2  *  *V  Sln  ^  ’ 

2(trr  ' 

cos  2qp  =  -  -  - - ^  , 

{ !|(  ffr =Crz  >  2+Crr  z  ^ 2 


where 


and 


-  13  ~ 


The  analogous  plastic  strain-rate  components  in  the  r  and  z  directions  are 


1  1 

€r  =  2  ^G1  +  €2^  +  2  ^€1  “  €2^  C°S  2<P 
€z  =  2  ^€1  +  62^  "  2  ^€1  “  €2^  C°S  29 
€rz  =  2  “  €2^  Sin  2<3p  s 


where  cos  2cp 


2  (V€z} 


{f(e  -e  )2+e2 
'•4  r  z'  rzJ 


1  9 


and 


(2.2.8) 


sin  2^ 


rz 


{f(  e  =e  )2+e2  )2 
u4  r  z'  rzJ 


The  condition  of  isotropy  is  conveniently  expressed  by  the  relationships 


CT 

r 


(2.2.9) 


which  follows  from  ( 2„2cY )  a^d  (2c2o8)o  The  condition  of  incompressibility 
is  defined -by- 


€1  +  €2  +  €3  ~  er  +  €z  +  £$  ~  (2o2ol0) 

The  general  field  equations  for  quasi-static,  axially  symmetric 

arc 

plastic  flow  for  the  type  of  material  considered  there  for  given  by 
equations  (2„l„l),  (20l„2)s  and  Table  I„  Equations  (2<,2„9)  and  (2„2„10) 
must  also  be  satisfied  by  the  stress  and  plastic  strain-rate  components,. 


r -•*/*  -  x 


-  Ik  - 


CHAPTER  III 

ANALYSIS  OF  THE  FIELD  EQUATIONS 

The  field  equations  of  Chapter  II  are  specialized  into  each  of  the 
various  plastic  regimes  of  groups  I  to  IV  and  a  mathematical  analysis  of 
their  structure  obtained „ 


3* 1  Group  lz  Plastic  Regimes  B  and  E, 


The  regimes  B  and  E  of  Fig0  2  are  singular  and  semi-isotropic 

since  cx,  =  cr  d  a  „  B  differs  from  E  only  in  that  it  corresponds  to  a 
1  d  '  3 

higher  mean  value  of  the  principal  stresses  for  some  given  value  of  the 


circumferential  stress  =  cr^ 


(a)  Stress  Fields„ 

The  yield  condition  for  B  (Table  l)  is  <r 1  =  a2  =  °3  +  Hence^ 

tj  *  a  -2k  since  by  (2o2<>7)s  <J  -  cr  and  a  =  0.  Therefore, 

6/  r  r  z  r  z 


dcr  2k 

+  ~~  s  0  by  (2.1.1)  and 


<j  =  o"  =  2k  In  "  »  A  >  0  s> 

r  z  r 

<T0  =  2k  [  in  (j)  -1]  . 


(3-1.1) 


This  stress  field  has  a  singularity  at  r  =  0  and,  consequently,  must  be 
defined  only  for  r  >  0  to  avoid  infinities  in  the  stresses. 


The  stress  field  for  E  can  be  derived  in  a  similar  manner;  the 
only  component  being  different  is  cr 

<rg  =  2k  [  in  (~)  +  1  ]  . 

1  - 


9 


where 


, 
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(b)  Velocity  Fields, 

From  Table  I,  the  plastic  strain-rate  components  for  B  are 


\  >  o. 


€2  =  —  0  and  =  -A^  -A^  <  o 


Since  €  =  s  “  by  (2,1,3)»  the  radial  velocity  u  must  be  zero  or 

negative  for  r  >  0„  There  being  no  shearing  stress  in  the  r  -  z  plane, 
the  principal  axes  of  stress  are  not  uniquely  determined  and  the  isotropy 

i 

condition  (2,2,9)  cannot  be  used  as  a  relationship  for  determining  the 


principal  strain-rates 0  Consequently,  the  velocity  field  is  not  deter- 

I  .  1 

minate,  being  governed  only  by  the  single  equation  of  incompressibility 
(2,2. 10),  The  velocity  components  u  and  w  must,  however,  be  contin¬ 
uous  functions  of  r  and  z,  This  follows  as  discontinuities  in 
tangential  velocities  can  occur  only  across  a  surface  on  which  the  shearing 
stress  is  k. 


The  velocity  field  at  E  is  similarly  indeterminate <,  The 
component  u,  however,  is  either  zero  or  positive  as  follows  from  Table  I, 


3„2  Group  Hi  Plastic  Regimes  AB  and  EF, 

The  plastic  regimes  AB  and  EF  are  regular.  Because  determination 
of  the  velocity  fields  and  stress  fields  for  these  regimes  is  similar, 
attention  is  confined  only  to  AB,  The  outstanding  difference  between  these 
two  regular  regimes  is  in  the  type  of  plastic  flow  possible  in  the  axial 
planes. 


From  Table  I,  the  principal  strain-rate  components  for  AB  are 

€ i  =  A i  >  0,  =  0,  and  =  -A^  <  0  , 

Hence,  e  *  =  -  <  0  and  this  implies  that  u  <  0  for  all  r  >  0,  The 

’  3  9  r  -  - 
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principal  strain-rate  components  for  EF  are 


€,  -  0,  €  s  -  A  <  0a  and  €  -  A„  >  0 

1  d  d  3  2““ 

Hence,  -  €0  “  "  >  0  and  this  implies  that  u  >  0  for  all  r  >  0.  The 
regime  AB  corresponds  to  states  of  stress  where  "waisting”  occurs  in  the 
axial  planes 9  in  contrast  to  the  regime  EF  where  "barrelling”  occurs  in 
the  axial  planes „ 


(a)  Velocity  Field  for  AB„ 

The  condition  that  -  0  requires  that 


€ 

r 


+  € 

Z 


€  )2  +  he2  P 

z ’  rz  J 


by  (2„2o6)»  Using  (2„ 1„3) »  this  can  be  written  as 


/  6u  6w  N\2  /  6w  f  6u  6w  \2 


Also  by  (2olo3)s  the  incompressibility  equation  (2.2ol0) 


(3»2ol) 

can  be  written  as 


6u  6w 

7Z  +  ^ 


+  H  =  o 

r 


(3-2.2) 


From  the  yield  condition  for  AB  (Table  X),  (T.  =  cr  +  2k,  <x  >  <r  >  <j  s 

*1  -  a2  3 

the  further  condition  that  -  ~  r"  <  k  shows  that  the  maximum  shearing 
stress  in  the  r  =  z  plane  must  be  less  than  k„  Consequently,  velocity 
components  u  and  w  are  to  be  continuous  functions  of  r  and  zc  These 
velocity  components  are  determined  by  (3.2„l)  and  (3*2<>2)  with  u  having 
the  additional  requirement  of  being  non-negative  for  all  r  >  0*  The 
velocity  field  for  AB  is  kinematically  determinate  in  the  sense  that  there 
is  available  as  many  equations  involving  the  velocity  components  as  there 
are  unknown  velocity  components „ 


The  equation  of  incompressibility  (3o2»2)  can  be  written 


-  IT 


(3-2.3) 


If  u  and  w  are  derivable  from  a  velocity  function  V(rsz)  such  that 


u  - 


1 

r  ciz 


and  w 


1  Sv 

r  Sr 


(3«2»4) 


then^ clearly,  V  satisfies  (3.2.3).  ¥  must  also  satisfy  (3. 2.1).  Sub¬ 

stitution  of  (5„2o^)  into  (3.2.I)  gives 


1  /  dv  V  /  1  dv  2  £%  V 

7  V  s  ;  =  -  72  sr 4  7  ss:  )  + 

which  can  be  simplified  to 


1  1  dv  1  V 

+  ’ 


f  1  dv  V  .  d%  /  ci^V  1  dv  \  ,  _  _  _ . 

^  57s  '  ‘  7  57 +1*S7S7(v§757-7S7j=0*  (5-2.5) 

a  non-linear  second-order  partial  differential  equation  in  V.  If  u9 
w  and  V  are  prescribed  as  initial  conditions  on  some  open  curve  F 
in  the  r  =  z  plane,  then  these  conditions  together  with  (3.2.5)  constitute 
a  general  Cauchy  problem  in  the  theory  of  partial  differential  equations 
(Petrovsky  [1961]).  The  classification  of  (3<>2.5)  and  the  possible 
extension  of  V  for  points  in  the  immediate,  neighbourhood  of  V  by  a 
Taylor  series  is  now  considered. 


The  following  relationships  must  be  satisfied  along  F  *  namely? 


d 

d 


dr  + 

■  dr  + 


d*V  , 
Sz2  dZ  • 


(3-2.6) 


Together  with  (3,, 2. 5)  these  relationships  are  sufficient  for  the  determin¬ 
ation  of  the  second-order  derivatives  of  V  provided  they  exist.  Similarly 
for  the  second-order  differential  coefficients  along  T,  there  are  the 


relationships 


ano 
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(3.2.7) 


In  addition^  there  is  the  relationship  obtained  by  differentiating  (3*2*5) 
partially  with  respect  to  r  yielding 
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(5°2o8) 

-  0 


a  quasi-linear  partial  differential  equation  of  the  third-order*  Equations 
(3*2*7)  and  (3»2*8)  are  sufficient  for  the  determination  of  the  third-order 
partial  derivatives  of  V  provided  the  coefficient  determinant  is  non-zero* 
Curves  for  which  this  determinant  is  zero  are  the  characteristics  of  the 
system  (3*2*5)  and  an  extension  of  V  is  not  possible  from  them  as  a  Taylor 
series  for  V  cannot  be  formed  (Schiffer  [i960]).  The  characteristics  of 
(3.2.5)  are  found  by  the  determinantal  equation 
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which  upon  expansion  becomes 
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(3.2.9) 
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_  (  5%  3% 

»r>z)  -  S3  -  gp 


Letting 
equation  (3.2.9)  becomes 


1 
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end  B 


*  !  ^ 
r,z)  =  U 


B 

dz2  -  —  dz  dr  -  dr2  =  0  . 


Now  (3.2.5)  can  be  written  as 


Substitution  of  (3.2. ll)  into  (3.2. 10)  yields 


dr  dz  -  dr2  =  0 


which  factors  into 


[• 
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Hence, either 


dz 


„  3% 

2  5757 


dr  "  3%  3^  1  3V 
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or 


2  £>V 
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(3.2.10) 


(3.2.11) 


(3.2.12) 


By  use  of  (3.2.5),  equation  (3.2.12)  can  be  written  as 
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The  characteristics  are 


thus  defined  by  the  equations 


and 
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(3.2.13a) 

(3.2.13b) 


The  partial  differential  equation  (3<>2o5)  is  therefore  hyperbolic.  Also 
the  two  families  of  characteristics  (joSdJa)  and  (3<,2„13b)  form  an  orthogonal 
net  since  the  product  of  their  respective  slopes  is  -  1  . 

From  Table  Is  the  plastic  strain-rates  for  AB  are 


x  =  \  >  0  ,  €2  =  0  »  €  =  ”A1  <  0 


Hence  by  (2o2c8) 
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by  (3.2.U)  , 
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by  (3.2.13a) 


Similarly^  since 
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dz  by  (302013b). 
dr 


But 
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Therefore 


=  tan  cp 
dr  T 


o 


(3.2. 1U) 


(3.2.15) 


From  (3.2.11+)  and  (3.2.15)>  it  is  seen  that  the  characteristic  directions 
are  the  same  as  those  of  the  principal  strain-rates  in  the  r  -  z  plane. 


(b)  Stress  Field  for  AB „ 

Structure  of  the  stress  field  for  plastic  regime  AB  is  determined 
from  the  equations  of  equilibrium  (2.1.1)  and  (2.1.2),  the  yield  condition 
from  Table  I,  and  the  isotropy  condition  (2.2.9). 


Consider  $(r,z)  as  a  stress  function  such  that 
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(3.2.16) 


$(r,z)  satisfies  the  equilibrium  equation  (2.1.2);;  namely., 
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“87“  +  ST 
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=  0 


The  isotropy  condition  (2.2.9)  can  be  written 


a  -  a 
r  z 


rz 


€  -  € 
r  z 


rz 


=  2  cot  2cp  by  (2.2.8) 


Therefore 


a  -  a  =  2cr  cot  2<P 
r  z  rz 

°r  =  r  '{  *2  COt  2(p)  §|  ‘  If  } 


(5.2.17) 


in  terms  of  the  stress  function  $(r,z).  The  yield  condition  to  be 
satisfied  is  <7^  -  O’  -  2k  or,  equivalent  ly,  by  use  of  (2.2.15)  that 

“  (cr  -  cr  )  +  {  f  (cr  -  cr  )2  +  a2  -  cr  =  2k  . 

2  r  z  ^  4  r  z’  rz  J  6 

Expressed  in  terms  of  $(r,z),  this  yield  condition  becomes 


1  8$ 


^  (  cot  2qp  +  esc  2cp  )  -  cr^  =  2k  . 


Hence 


1  8$ 
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9  r  oz 


(3.2. 18) 


The  equilibrium  equation  (2.1.1)  expressed  in  terms  of  $(r,z)  now  becomes 


82$  82$  cot  cp  8$  , 

-  5?  -  2  cot  2<p^  +  — P*  ^  -  2k  =  0  , 


(5.2.19) 
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a  linear  second-order  partial  differential  equation  in  <£>.  The  characteristic 
directions  are  determined  using  (3.2,19)  and  equations  (3.2,20);  namely 
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which  expresses  the  variation  of  the  first  derivatives  of  0  along  any 
curve  T  on  the  rsz  plane.  Equating  to  zero  the  coefficient  determinant 
of  (3,2,19)  and  (3,2,20)*,  ie. 


t* 

1 

-2  cot  2cp  ■ 
4 

-1 

»  dr  >' 

dz  ; 

0 

ii 

0 

0 

dr 

dz 

yields  upon  expansion 

dz2  +  2  cot  2cp  dr  dz  -  dr2  =  0  . 


dz 

Thus  — —  =  tan  cp  * 

dr  i 

dz 

and  *7“  =  -  cot  cp  are  the  characteristic  directions, 

dr 

The  above  analysis  shows  that  the  partial  differential  equation 
(3,2.19)  is  hyperbolic  and  that  the  characteristics  form  an  orthogonal  net 
coincident  with  the  characteristics  of  the  partial  differential  equation 
(3°2,5)  the  velocity  function  V(r,z), 


In  conclusionp  therefore,  the  plastic  regime  AB  (similarly  for 
EF)  has  a  velocity  field  which  is  kinematically  determinate  and  hyperbolic 
and  a  stress  field  which  is  statically  determinate  and  hyperbolic.  The 
characteristics  of  both  fields  are  coincident  and  have  the  same  directions 
as  the  principal  strain-rates  in  the  r,z  plane. 


■ 
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3°3  Group  Ills  Plastic  Regimes  A  and  F„ 

The  plastic  regimes  A  and  F  in  Figure  2  are  singular  and 
characterized  by  the  equality  of  the  circumferential  principal  stress  with 
one  of  the  two  principal  stresses  in  the  meridian  plane.  This  is  the  "full 
plasticity  (vollplastisch)"  hypothesis  proposed  by  Haar  and  von  Kerman  in 
1909. 

(a)  Stress  Field  for  F. 

Consider  regime  F  for  definiteness.  There  are  two  equilibrium 

equations  (2,l,l)  and  (2,1,2)  and  two  yield  conditions  from  Table  I  available 

for  the  determination  of  the  four  stress  components  cr  ,  a  ,  O'  and  or 

r  0  z  r  z 

In  this  sense^the  stress  field  is  statically  determinate. 


From  the  yield  conditions  or^~ 
it  follows  that 


2k  and 


2k  where 


by  the  use  of  (2,2.5). 


Hence  the  maximum  shearing  stress  in  the  meridian  plane  is  k.  The  yield 
conditions  impose  two  independent  conditions  on  the  four  stress  components 
and3 consequent ly; the  possible  yield  states  may  be  represented  by  two  inde¬ 
pendent  parameters.  These  parameters  are  taken  as  the  pressure 
p(r.z)  =  -  |  (0^  +  Og) 


(3.3.1) 


and  the  angle  q>(rS)z)  which  specifies  the  orientation  of  the  principal 
axes  in  the  meridian  plane.  The  lines  of  maximum  shearing  stress  (slip 
lines)  form  an  orthogonal  system  of  curvilinear  coordinates  at  a  point  in 
the  meridian  plane.  By  convention,  these  lines  are  called  a-  and  £  -lines 
and  the  direction  of  the  algebraically  greater  principal  stress  0^  is 


' 

/ 


obtained  by  an  anticlockwise  rotation  of  an  angle  at  A  from  the  a  -  lineo 
The  inclination  of  the  a  «■  line  at  any  point  with  the  r-axis  is  <p(rsz)„ 

a»  (B  -  lines  are  acted  upon 

-  P  »  ' 

and 

crr  =  =  p  -  k  sia  2(P) 


Since  surface  elements  perpendicular  to  the 
by  a  shearing  stress  k  and  a  normal  stress 


a 


a 


3 


<j  -  =  p  +  k  siu  2®, 
z  ’ 


a  -  k  cos  2cp  , 

r  z  ’ 


(3.5.2) 


%  =  -  p  + k  • 

Substitution  of  (3.3.2)  into  the  equations  of  equilibrium  (2.1.1)  and 
(2.I.2)  results  in 


Op  0Q) 

+  2k  cos  2cp  ^  +  2k  sin  2qp 


(1  +  sin  2cp) 


d  d  S  1c 

and  +  2k  sin  2<p  ^  2k  cos  2q>  “  cos  2q> 


(3.3.3) 


The  characteristics  of  this  system  of  quasi-linear  simultaneous  equations 
are  determined  using  (3°3°3)  together  with  the  relations 


(3.3A) 


by  equating  to  zero  the  coefficient  determinant  of  Cramer's  rule  for  the 
dp  dp  dt 


solution  of  ^  and  ^  (Hildebrand  [195^]) -  Specifically, 


all  - 


'•  ’V  -  v  ‘  ■*  1 


=*  26  •= 

1  0  2k  cos  2cp 
0  1  2k  sib  Sep 

dr  dz  0 

0  0  dr 

yields  upon  expansion  the  equation 

(sin  2cp)  dr2  -  2  cos  2cp  dr  dz  ~  (sin  2cp)  dz2  =  0  . 

Therefore 

dz 

^  s  tan  cp 

(3.3.5) 

dz 

or  dr  ~  °  COt  ^  0 

Accordingly,  the  characteristics  (3<,3<>5)  are  teal  and  distinct  and  the 
stress  equations  (3°3°3)  arft  hyperbolic0  The  characteristics  are  orthogonal 
and^moreover^ coincide  with  the  lines  of  maximum  shearing  stress 
(a  and  (3  ■=»  lines) „ 


2k  sin  2cp 
>2k  cos  2cp 
0 

dz 


=  0 


Requirements  to  be  satisfied  along  the  characteristics  (3<>3°5) 

6  6  6  6 

to  insure  solutions  for  ^  and  ^  are  obtained  by  equating  to 

zero  the  numerator  determinant  of  Cramer's  rule0  Specifically,  this 
determinantal  equation 


1  0  2k  cos  2q) 

0  1  2k  sin  2q> 

dr  dz  0 


( 1+sin  2q>) 
“  (cos  2cp) 
dp 


0  0  dr 


dq> 


=  0 


yields  upon  expansion 


/ 
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d  z  1c 

dp  +  JJT  (2k  sin  2cp  dcp  -  ~  cos  2cp  dr) 


4  2k  cos  2cp  dcp  +  “  (l  +  sin  2cp)  dr  =  0 


But  for  an  a  °  line 


dz 

dr 


tan  cp9  dr  =  cos  cp  ds  and  dz 


sin  cp  ds  , 

CL 


(3.3.7) 


where  ds  is  the  differential  arc-length.  Thus  for  an  a  -  line. there 

CL  f 

is  the  requirement  that 


dp  4  2k  dcp  +  k(sin  cp  +  cos  cp) 


a 


0 


(3.3.8) 


Similarly  for  a  (3  -  line. 


dz 

dr 


cot  cp j>  dr  = 


sin  cp  ds^  and  dz  =  cos  cp  ds^  ,  (3.3.9) 


where  ds  is  the  differential  arc-length.  Accordingly  for  a  3  -  line 
P 

there  is  the  requirement  that 


ds 


dp  =  2k  dq>  -  k  (coscp  +  sin  cp)  =  0 


(3.3.10) 


(b)  Velocity  Field  for  Fc 

For  plastic  regime  F,  the  velocity  components  u  and  w  are 
determined  from  the  isotropy  condition  (2.2.9)  and  the  incompressibility 
condition  (2.2.10)  or  (3.2.2)  as  it  is  more  conveniently  written  for 
purposes  of  this  analysis.  The  isotropy  condition  to  be  satisfied  is 


as  follows  from  (3.3.2)  and  (2.2o9)<>  There  are?  therefore^  two  equations 
(3.2.2)  and  (3.3.II)  available  for  the  determination  of  u  and  w; 


cot  2cp 


(3.3.11) 


namely, 


' 


*alas?9iab  arfj 
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*.  rw  du  ^  _  8w 

cot  351  57  +  51  +  57  -  cot  2cf  57 


«  0 


and 


8u  8w  u 

+ 

r 


(3.3.12) 


57 + 57 + ~ = ° 


Also? there  exists  the  restrictions  from  Table  I  that  <  0 


u 


and  e,  =  K  >  0  .  Since  e _  =  e _ 

3  2“  3  e  r 


u  >  0  for  r  >  0.  Also 


<  0  requires  that 

^  X.T  \2 


/  8u  8w 

(57-57 


/  8u  8w  N2  1 2  u 

+  \  si +  ;  1  *  "  7 


as  follows  from  (2.103)j  (2»2o6)  and  (302„2) 


The  characteristics  of  the  velocity  field  are  determined  from 
(3o3„12)  and  the  additional  relationships  that 


and 


,  8u  ,  8u  j 

du  =  57  dr  +  57  dz 

.  8w  ,  8w  , 

dw  =  ^  dr  +  ^  dz 


(3°3° 13) 


for  any  open  curve  F  on  the  rsz  plane  along  which  u  and  w  are  known, 
These  characteristics  are  defined  by  the  determinantal  equation 


cot  2cp 


dr 


0 


0 


dz  0 


0  dr 


which  upon  expansion  yields 


cot  2cp 
1 
0 
dz 


dz2  +  2  cot  2cp  dr  dz  -  dr2  =  0 


=  0 


Therefore  “  =  tan  cp 
dr 


and 


dz 

dr 


COt  Cp 


(3.3.1U) 
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Hence  the  characteristics  of  the  velocity  field  are  orthogonal  and  coincide 
with  the  characteristics  of  the  stress  fieldo 


The  relationships  along  these  characteristics  are  obtained  from 


the  determinantal  equation 


This  reduces  to 


u 

r 

du 


1 

0 

dz 


1 

0 

0 


dw  0  dr 


cot  2 © 

1 

0 

dz 


=  0 


”  dz2  +  dw  dz  +  du  dr  +  “  cot  2©  dz  dr  =  0  . 
r  r 


(3.3.15) 


By  use  of  (3.3.T)»  the  relationship  along  an  a  *  characteristic  is 


ds 


cos  ©  du  +  sin  ©  dw  +  “  —Sk  _  q 


r  z 


(3.3.1 6) 


Similarly  by  (3°3°9)»  the  relationship  along  a  3  -  characteristic  is 


sin  ©  du  -  cos  ©  dw 


u 

r 


ds 

z 


*  0 


(3.3.17) 


These  relationships  (3<,3„l6)  and  (3°3°17)  are  conveniently  written  in 
terms  of  the  velocity  resolutes  U  and  Wc  If  U  is  the  velocity  component 
along  an  a  -  characteristic  and  W  is  the  velocity  component  along  a 
3  -  characteristics,  then 


U  s  u  cos  ©  +  w  sin  ©  9 

and  W  =  -  us  in©  +  w  cos  (p  , 


(3.3.18) 


From  (3o3o18)s  it  follows  that 


,  p  II  f  * 

, 


“  30  =* 


dU  ss  cos  cp  du  +  sin  q>  dw  +  Wdq> 
and  dW  =  sin  cp  du  +  cos  qp  dw  -  Udcp  „ 


(3.5.19) 


The  characteristic  relations  (3.3. 16)  and  (3.3. 17)  now  reduce  to 

dU  -  Wdcp  +  ~  ds^  =  0  on  an  a  -  line,  (3.3,20) 

and  dW  +  Udcp  +  ^  d$^  =  0  on  an  p  -  line.  (3o302l) 

In  conclusions  the  plastic  regime  F  is  an  application  of  the 
Haar  and  von  Karman  hypothesis.  The  stress  field  is  statically  determinate 
and  hyperbolic  with  characteristics  coincident  with  the  lines  of  maximum 
shearing  stress.  The  velocity  field  is  kinematically  determinate  and 
likewise  hyperbolic  with  characteristics  identical  with  those  of  the  stress 
field.  These  statements  apply  to  the  plastic  regime  A  also. 


3<>4  Group  IV °  Plastic  Regime  AF. 

(a)  Velocity  Field  for  AF. 

The  plastic  strain-rate  components  for  this  regular  regime  are 


€ ,  -  A  >  0  s 
1  3 


€  =  -  A  <  0  and  e_  =  0 

^  3  ”  3 


u 


as  found  in  Table  I.  Since  “  =  0,  there,  is  no  displacement  or 


Su 


velocity  component  u  in  the  radial  direction.  Moreover^ since  ^ 

=  0.  From  the  incompressibility  equation  (3°3»2)9  it  follows  that 

£  =  s  0  .  Hence  €  -  -  is  the  only  non-vanishing  plastic  strain- 

z  Oz  rz  ox* 

rate  component  and  for  this  w  must  be  a  function  of  r  only.  With  € 
being  the  only  non-zero  component  it  means  that  the  r  and  z  directions 
are  the  directions  of  the  maximum  shearing  strain-rate. 


* 
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(b)  Stress  Field  for  AF0 

For  the  determination  of  the  stresses B  there  is  available  the 
yield  criterion  cr,  =  +  2k  where  cr.  >  cr  >  <t0,  the  isotropy  condition 

12  132 

(2o2o9)  and  the  equations  of  equilibrium  {2.1.1)  and  (2.L2),  The  yield 


condition  implies  that 


s  k  0  Thus  the  maximum  shearing  stress  in 


the.  meridian  plane  is  k.  The  isotropy  condition  (2.2.9)  becomes 


±  k 

c 

rz 


rz 


rz 


and 


(3.4.1) 

(3.4.2) 


From  (3o4o l)5 choose  cr^  -  +  k  for  definiteness.  Equation  of  equilibrium 
(2olo2)  becomes 


5T 


z  +  * .  0 

r 


Hence 


cr  -  cr  - 
z  r 


“  z  +  f(r)  » 


(3°bo3) 

(3ok,b) 


where  f(r)  is  an  arbitrary  function  of  r0  By  substitution  of  cr  into 


(2.1.1),  there  arises  the  equation 


Z  + 


f°(r)  -  “2  z  + 

'  /  r  -  r 


ikl.fa  . 


=  0 


Hence 


or 


<r0  =  r  f'(r)  +  f(r) 


°e  ■  3?  [r  f(r)1 


(3-4.5) 


The  function  f(r)  must  be  such  that  cr&  is  an  intermediate  principal 
stress.  From  (20205)>this  requires  that 


|  (Tr  +  az)  +  k  >  a0  >  g  (<rr  +  »z)  *  k 


or?  equivalently^  that 


' 
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-  r  2  +  k  >  r  I?  f<r)  >  -  7  1  ‘  k  •  (3*^*6) 

Thus  for  plastic  regime  AF?  the  directions  of  the  principal 
strain-rates  and  the  principal  stresses  are  fixed*  The  stress  field  is 
statically  determinate  and  the  velocity  field  is  kinematically  determinate* 
The  fact  that  the  radial  velocity  component  is  zero  renders  the  application 
of  this  regime  to  limited  and  rare  cases* 


' 
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CHAPTER  IV 

INCIPIENT  PLASTIC  FLOW  IN  A  S EMI- INFINITE  REGION  OF  RIGID-PLASTIC 
NON -HARDENING  ISOTROPIC  MATERIAL  DUE  TO  A  LOAD  APPLIED  BY  A  SMOOTH 
CONICAL  XNDENTER  FITTED  INTO  A  CONICAL  CAVITY 

1 

The  problem  of  finding  an  equilibrium  plastic  stress  field  under 
axially  symmetric  conditions  for  indentation  by  a  cone  is  considered .  The 
yield  point  load,  which  is  the  load  at  which  the  material  first  deforms,  is 
derived  using  the  incomplete  stress  field 0  No  bounds  are  made  on  the  actual 
yield  point  load  since  the  velocity  field  and  the  extended  stress  field  into 
the  rigid  region  are  not  determined 0 

U  o 1  Preliminary, 

The  adopted  procedure  used  In  this  problem  is  similar  to  that 
used  by  Berezeneov  [1955]*  Shield  [1955]*  and  Cox,  Eason  and  Hopkins  [1961], 
Each  investigator  assumed  a  priori  a  plastic  regime  applicable  to  their 
respective  problem,,  Except  for  Berezeneov,  they  verified  that  the  derived 
plastic  stress  field  was  the  correct  one  by  determining  a  compatible  kin¬ 
ematically  admissible  velocity  field.  The  validity  of  this  procedure  has 
been  established  by  Hill  [1951]  and  Bishop  [ 19533 » 

Hill  [19513  showed  by  virtual  work  and  the  maximum  work  principle 
for  an  element  that  wherever  deformation  is  occurring  in  a  rigid-plastic 
non-hardening  materials  the  state  of  stress  is  unique  in  the  deforming 
regions  of  solutions.  In  the  rigid  zone  common  to  all  the  solutions,,  the 
stress  need  not  be  unique.  Similar  conclusions  could  not  be  drawn  for  the 
velocity  field,  in  fact.  It  need  not  be  unique  anywhere.  The  velocity  field 
must  only  be  compatible  with  the  stress  field.  Hill  [1950c]  has  also  shown 
that  a  correct  velocity  field  is  obtained  in  the  plastic  zone  if  upon 
examination  of  the  stress  increment  occurring  during  an  increment  of 


I 
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distortion;,  the  elements  assumed  deformed  remain  plastically  stressed,. 

Actual  velocity  fields  occurring  in  real  metals  depend  upon  factors  such  as 
elastic  effects  and  work-hardening „  These  factors  are  neglected  in  ealeul- 
ations  for  the  yield  Xo,ad  on  rigid-plastic s  non-hardening  metals „ 

Bishop  [19533  discussed  a  working  procedure  for  the  solutions  of 
two-dimensional  kinematically  determinate  plastic  problems 0  The  inverse 
of  Bishop's  procedure*,  as  used  by  Shield  [19553 9  is  more  applicable  to 
problems  that  are  statically  determinate,.  Under  this  scheme,,  one  proceeds 
by  ignoring  the  development  of  the  plastic  zone  and  derives  the  unique 
plastic  stress  field  and  an  associated  distortion  mode  that  satisfies  the 
boundary  conditions  at  yield  only,.  This  forms  an  ’’incomplete"  solution  in 
Bishop's  terminology,,  For  a  ’’complete"  solution  a  stress  field  is  required 
in  the  rigid  zone»  The  stress  distributions  must  satisfy  the  equilibrium 
equations  and  must  not  violate  the  yield  condition,,  According  to  Hill 
[ I95I ] 9  if  the  yield  condition  is  locally  exceeded;,  such  a  solution  is 
applicable  for  a  metal  which  is  at  least  hardened  there  by  that  amount . 

Extremum  principles  established  by  Harkov  [19^-7 3 »  Hill  [i95^d, 

I95I ] 9  Prager  [I95U]  and  limit  analysis  theorems  (Drucker,  Greenberger 
and  Prager  [1952])  can  be  used  on  rigid-plasti^ non-hardening  materials 
to  find  bounds  for  the  yield  point  load,.  Hill  [19513  Has  shown  that  any 
distribution  of  stress  satisfying  the  equilibrium  equations,  the  stress 
boundary  conditions  and  nowhere  violating  the  yield  criterion  gives  rise 
to  external  boundary  loads  which  are  not  greater  than  the  actual  yield  point 
loado  The  partial  stress  field  obtained  in  an  incomplete  solution  does  not 
meet  the  requirements  of  this  lower  bound  theorem  , consequently, cannot  be 
used  to  establish  a  lower  bound,,  However,  if  it  can  be  demonstrated  that  an 

field  (not  necessarily  unique)  that  satisfies  the 


equilibrium  stress 


' 
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boundary  conditions  and  does  not  exceed  the  yield  point  exists  in  the  rigid 
region9  then  this  complete  stress  field  can  be  used  to  establish  a  lower 
bound o  The  extremum  principles  also  show  that  the  internal  rate  of  plastic 
work  calculated  from  any  velocity  field  or  mode  of  deformation  compatible 
with  any  plastic  stress  field  and  velocity  boundary  condition  is  not  less 
than  the  rate  of  working  of  the  actual  boundary  stresses „  Any  postulated 
velocity  field^ therefore^ gives  an  upper  bound  to  the  actual  yield  point 
loado 


Shield  [19553  used  the  above  procedure  to  solve  the  Ishlinskii 
problem  referred  to  in  Chapter  Ic  By  assuming  that  the  state  of  stress  at 
yield  is  represented  by  the  Haar-Karman  plastic  regime  F  of  Group  III, 
Chapter  III,  Shield  derived  a  complete  solution  and  verified  it  as  the 
actual  oneD  This  was  attained  by  showing  that  the  upper  and  lower  bounds 
for  the  yield  point  load  were  equal 0 

4,2  Formulation  of  the  Problem, 

With  reference  to  Figure  3,  for  a  cylindrical  polar  coordinate 
system  (r,  0,  z)  with  the  z°axes  ehoosen  as  vertical,  let  the  region 
z  >  0  be  occupied  by  a  rigid-plastic^non-hardening,  isotropic  material 
except  for  the  region  where  a  right  circular,  conical  cavity  with  45° 
semi-angle  exists  such  that  its  base  of  radius  R  is  on  the  surface 
z  =  0.  The  axis  of  the  cavity  is  perpendicular  to  the  surface  of  the 
material  and  the  center  of  the  base  is  ehoosen  as  the  origin  of  the 
coordinate  systemD  A  smooth,  rigid,  right  circular,  conical  punch  is 
inserted  into  the  cavity  fitting  it  fully„  The  punch  is  then  centrally 
loaded  until  incipient  plastic  flow  of  the  material  occurs.  The  arrange- 
ment  is  clearly  one  of  axial  symmetry. 
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For  axial  symmetry  [Chapter  II  (2„l)]s  the  circumferential  stress 
cr  must  be  a  principal  stressc  The  surface  of  the  material  for  which  r  >  R 

v7 

is  stress  free,  Consequently  or  =  or  =  0  on  this  surface0  This  condition 

z  rz 

requires  that  and  or^  be  principal  stresses  for  surface  elements  and 

that  in  any  meridian  plane  the  lines  of  maximum  shearing  stress  meet  the 
free  boundary  at  1+5°  angles  0 

The  state  of  stress  at  incipient  plastic  flow  is  assumed  to  be 
represented  by  the  Haar-Karman  plastic  regime  F  of  Group  III  (Chapter  IXX)„ 
Consequently^ the  equations  governing  the  stress  state  in  the  plastic  region 
are  hyperbolic  and  the  a  °  and  jB  -  characteristics  are  given  by  (3°3°5)° 

dz  f  tan  q>  on  an  a  -  line  , 

dr  ~  [  -cot  cp  on  a  p  -  line  0 

The  relationships  along  the  characteristics  are  given  by  (3°3°6)  anc*  (3°3°10)° 

ds 

dp  +  2k  dep  +  k(sin  cp  +  cos  q>)  =  0  on  an  a  -  line  , 

ds 

and  dp  -  2k  dcp  -  k(sin  cp  +  cos  cp)  =  0  on  a  p  -  line  . 

In  Chapter  III,  the  characteristics  of  this  plastic  regime  were  shown  to 
coincide  with  the  lines  of  maximum  shearing  stress „  Hence  they  must  meet 
the  stress  free  boundary  z  =  0,  r  >  R  at  45  angles 0 

Since  the  development  of  the  plastic  regime  is  ignored,  let  the 
region  ABE  of  Figc  3  represent  one-half  of  the  plastic  stress  region  in  a 
meridional  plane  at  yield,  the  remaining  half  of  the  field  follows  by 
symmetry,, 

4o3  Plastic  Stress  Field  ABC. 

From  hyperbolic  characteristic  theory  [Schiffer,  i960],  the  state 
of  stress  can  be  determined  in  the  region  ABC  (ligo  4)  from  the  initial 
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conditions  on  the  non-character istic  line  AB,  Line  AC  is  an  a  -  char- 
acteristic  through  singular  end  point  A  and  BC  is  a  j3  -  characteristic 
through  end-point  Bc  The  position  of  B  (Fig,  3)  on  the  stress-free 
boundary  is  not  known  until  it  can  be  determined  where  B  must  be  such 
that  an  extended  0  -  line  through  B  will  pass  through  the  E,  the  apex 
of  the  conical  cavity,  The  solution  of  the  stress  field  in  ABC  is  a 
Cauchy  problem  in  partial  differential  equation  theory  and  corresponds  to 
the  second  boundary-value  problem  [Hill,  I950e]  in  plane  strain  plasticity,, 

Shield  [19553  has  derived  general  equations  for  the  characteristic 
fields  generated  by  straight-line,  stress-free  boundaries  which  are  inclined 
at  an  angle  7  to  the  r  -  axis„  In  particular,  plastic  stress  fields 
were  obtained  for  7  s=  0°,  30°  and  60°,  The  field  for  7=0°  is  immed¬ 
iately  applicable  to  the  present  problem  and  an  analysis  of  Shield's  work 
is  considered  here  for  completeness „ 

A  straight-line  boundary  in  the  r  -  z  plane  of  an  axially 
symmetric  body  arises  when  part  of  the  surface  of  the  body  is  a  portion  of 
a  circular  cone0  Figure  5  illustrates  the  arrangement  whereby  a  cylindrical 
polar  coordinate  system  is  choosen  such  that  its  origin  coincides  with  the 
virtual  apex  of  the  ccnical  free  surface„  The  material  bordering  a  straight- 
line,  stress-free  boundary  is  assumed  t©  be  in  a  plastic  state  of  stress 
represented  by  point  F  of  Fig0  2  in  principal  stress  space. 

Let  ip  s  tan  ™  ~  ^  J  >  (U.3.1) 

and  the  free  boundary  be  ijf  ~  7  ,  7  a  constant.  Slip  lines  meet  the  free 

1  o 

boundary  at  45  angles. 

The  field  generated  by  the  free  boundary  is  determined  on  the 
assumption  that  p  and  cp  are  functions  of  ip  only;  an  assumption 
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compatible  with  the  equilibrium  equations  since,  no  fundamental  length  is 
involved,.  From  condition  (l4-„3°l)s>  the  operators 


d  1  2  •  d 


51  ■  7  cos  * 


dip 


and 


5? 


1  .  ,  ,  d 

“  Sin  Ip  cos  ip 


dtp 


(4.3.2) 


are  derivable,,  When  (3°3°2)  are  substituted  into  equilibrium  equations 
(2.1.1)  and  (2.1.2)  and  by  use  of  the  operators  (4o3„2)9  there  results 
the  differential  equations; 

k  d^  +  ^  d^  jcos  2cp  -  cot  ip  sin  2Kpj-  -  -jl  +  sin  2cpj-  esc  ip  sec  ip  =  0,  (4. 3° 3) 
°  k  +  2  dtyv  '|cos  29  +  tan  ip  sin  2cpj»  +  cos  2ep  sec ^  s  0,  (k«3ok) 

Adding  (4„3°3)  and  (4.3»^)  yields 

2  ^  -^2  cos  2q)  +  sin  2cp  (tan  f  =>  cot  ip )|»  (4o3<>5) 

+  cos  2q?  sec2^  -  (1  +  sin  2cp)  sec  ip  cosec  ip  -  0, 


Let 


then 


r  s  2(9  -  f)  > 


2  M  „  dT 

dip  ~  dip  * 


(4*3°6) 


Substitution  of  (4<,3°6)  into  (4o3°5)  yields 


dr 

dip 


sin  F  +  3sin  F  +  cos  F  tan  ^+1  =  0 


C^3°7) 


From  (4.3.4), 


“  ^  =  2  4?  ieos  2cp  +  tan  ip  sin  29b  +  cos  2q)  sec ^tp 
k  dip  dip  {  J 


By  use  of  (403°6),  this  becomes 


l  d£ 
k  dip 


'dT 

'dT 

'dip 


+  2 


+  3 


){ 


eoiiCr  +  ip) 

cos  ip 

cos(F  4-  #) 

cos  ^ 


cos  2cp 


cos 


sin(F  +  ip)  sin  ^ 

cos2^ 


, bavlovni 


dr 

dtp 


cos  r  +  cos  F  «  sin  F  tan  ip  +  4^)=cos  F  tan2^/. 


But 


dF  ,  ,  1  cos  F  tan  ifr  „  x 

+  4^1=  -T-^r  -  - - r— »-  by  (4  =  3.7) 


dip 


sin  F 


sin  F 


Therefore 


I  M  dF 
k  dip  ~  dip 


Hence 


2M 

k 


cos  F  +  3  cos  F  °  sin  F  tan  ip  +  tan  ip 


ip 


jT  (3  cos  F  -  sin  T  tan  ip)  dtp  +  sin  T  -  In  cos  ^  +  A 
7 

(M°8) 


where  A  is  determined  from  the  condition  that  p  =  k  when  If  =  7 


To  plot  the  a  -  characteristics,,  let  R  be  the  distance  from 

a 

the  origin  of  a  point  on  an  a  -  line  and  R  the  value  of  R  on  the 

ao  a 

boundary  ip  s  7  „  With  reference  to  Figure  5a»  there  is  derived  the 
relationship 


dR 


a 


R  dip 
tx 


tan  [  |  -  (cp  -  ip)  ] 


=  cot  (cp  -  ip) . 


ip 


Therefore  in 


or 


—  =  f  cot  (cp  -  ip)  dtp 
■„n  ^ 

f  f  cot  I  ^  } ' 


7 

exp 


(4.5.9) 


Similar ly? with  reference  to  Figure  5b ,  if  R^  is  the  distance  from  the 
origin  of  a  point  on  a  (3  -  line  and  R^q  is  the  value  of  Rr  on  the 


0 


boundary  ip  =  7  »  then 


=  tan  [  °(cp  -  ] 


Therefore 


V* 

R  1 
B  =  exp 

V 


tan 


} 


(4.3.10) 


For  the  present  problem  in  conical  indentation,  7=0  and  the 

boundary  conditions  on  ijj  -  y  =  0  are  p(y)  -  k  and  cp  =  jr/4  .  Using 

numerical  integration;,  Shield  determined  F(^)  from  (4.3.7)  together  with 

the  initial  condition  F(y)  =s  jr/2.  Similarly  from  (4.3.8),  p(^)  was 

R  Rfl 

obtained  with  the  initial  condition  that  p(/)  =  k,  and  -  were 

ao  30 

then  determined  from  (4„3°9)  and  (4.3°10)  by  numerical  integration.  The 


results  are  listed  here  in  Table  II0  Inspection  of  (4.3.7)  reveals  that 


dF 

becomes  undefined  when 
atf/ 

F  =  0.  When 

this  occurs, 

=  f  .  This 

means  that 

the  a°  lines 

approach  asymptotically  a  limit 

line.  This  limit 

line  is  a 

straight  line 

f  - 

=  0.301  radians. 

p  the  last  entry 

of  tne  first. 

column  of 

Table  II. 

Table  II 

if 

( radian) 

<P 

p/k 

R 

a 

Rao 

V 

V 

0.00 

0.7854 

1.0000 

1.0000 

1.0000 

0.02 

Oo7654 

0.9992 

1.0210 

0.9810 

0.04 

0.7453 

0.9968 

1.0443 

0.9637 

0.06 

0.7250 

0.9927 

1.0702 

0.9481 

0.08 

0.70421 

0.9868 

1.0991 

0.9339 

0.10 

0.6833 

0.9791 

1.1314 

0.9211 

0.12 

0.6617 

0.9693 

1 . 1679 

0.9096 

0.14 

0.6393 

0.9572 

1.2094 

0.8992 

0.16 

0.6160 

0.9423 

1-2571 

0.8900 

0.18 

0.5914 

0.9243 

1.3126 

0.8817 

0.20 

0.5652 

0.9023 

1 . 3784 

0.8746 

0.22 

0.5367 

0.8751 

1.4584 

0.8684 

0.24 

0.5048 

0.8408 

1 0 5594 

O.8652 

0.26 

0.4678 

0.7954 

1.6948 

0.8590 

0.28 

0.4207 

0.7293 

1.9005 

0.8560 

<Z2  <Z> 

0.301 

e  c» 

0.301 

°  00 

00 

0.854 
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*  and  in  Table  II  were  used  to  draw 

K  K 

oo  (3o 

the  plastic  stress  field  in  region  ABC  (Fig.  k)  <>  The  drawing  was  done  with 


The  values  of  if/ 


a  Nestler  drafting  machine „  R  for  the  a  -  line  AC  was  chosen  as 

ao 

R  =  25  cm.*  the  maximum  radius  of  the  conical  cavity,,  R  was  then 

ao 

incremented  by  k  cm.  for  each  successive  a  -  line.  Similar  procedure  was 
followed  for  the  (3  =>  lines  where  R^q  for  the  degenerate  (3  ~  line  at  A 
was  likewise  R  =  25  cm..  Figure  3  illustrates  only  that  region*  ABC, 
of  the  stress  field  calculated  in  the  above  manner  that  is  an  actual  part 
of  the  true  stress  field  for  the  indentation  problem  under  analysis. 
Originally  the  plastic  stress  field  had  been  determined  to  a  larger  extent 
as  in  Figure  l+s  but  this  was  reduced  only  after  the  plastic  stress  field 
had  been  extended  into  the  remaining  region  and  the  specific  (3  =>  line 


through  E  could  be  determined, 


koh  Plastic  Stress  Field  ACD. 

The  calculation  of  the  plastic  stress  field  in  region  AC^D 
(Figure  k)  is  obtained  from  the  boundary  conditions  in  the  neighbourhood 
of  the  point  A  and  from  the  values  of  p  and  cp  at  the  points  of 
intersection  of  the  a  °  characteristic  AC^  by  the  tJj  -  lines  listed  in 
Table  II.  These  points  of  intersection  did  not  coincide  with  the  points 
of  intersection  of  the  a  -  line  AC  by  the  (3  -  lines  determined  in  ABC. 
Because  p  and  q>  are  calculated  values  at  the  former  points*  they  are 
used  for  the  extension  of  the  stress  field  into  AC^D.  A  is  a  singularity 
through  which  all  the  a  -  characteristics  in  A^D  pass.  The  angular  span 
of  this  fan  of  a  -  characteristics  is  k5° »  determined  by  the  condition 
that  the  a  -  characteristic  AD  must  meet  the  smooth  boundary  of  the 
conical  cavity  at  a  k5°  angle.  The  arrangement  is  analogous  to  the  first 
boundary  -  value  problem  [Hill*  1950  ]  in  plane  strain  plasticity. 


' 


Consider  3  -  characteristic  PQ  of  lengths  in  Figure  4.  The 
values  of  p  and  cp  along  this  curve  are  governed  by  (3.30IO);  namely , 


dp  -  2kdqp  -  k(sin  cp  +  cos 


dSn 

— i  - 


0 


Integration  of  this  equation  along  PQ  is  represented  by 

P~  QK  S 

p  Q  p  Q  p 

/  dp  “  2k  /  dcp  =  k  /  (sin  cp  +  cos  cp)  — ~  =  0 
P  Cp  o 


But  dr  = 


■=  sin  cp  dSp  for  3  »  characteristics.  Therefore 

p  Q  j 

Pq  -  Pp  "  2k  (cpQ  -  cpp)  +  k  J  (1  +  cot  cp)  =  0.  (4.4.1) 


If  PQ  is  allowed  to  become  indefinitely  small  by  approaching  A,  then 
r  R  „  But  Pq,  Pp,  q>Q  and  cp^  will  each  have  distinct  values 

denoted  as  PA’,  PA»  cp'A  and  cp^  respectively. 


Now 


lim 
r  ->  R 


x* 

J  ^  (l  +  cot  cp)  =  0  since  1  +  cot  cp  is  bounded 


for  all  values  of  cp}(0.301  <  cp  <  it/2)*  in  the  re8lon  ACD°  In  the  limiting 
form5(4.4.  l)  becomes 


PA  2kCpA  =  PA 


2kcp . 


or 


V 


k  +  “V  •  ^a’ 


<PA  -9, 


(4.4.2) 


Since  Q  can  arbitrarily  be  any  point  Q.  on  3  -  characteristic  in 
AC, D  ,  (4. It. 2)  may  be  written 

(4.4.3) 


From 


A  i  =  fA  +  2^pA"  ,  ^PA°  -  cpA" 

k  k  i  i  i 

P 


Table  II,  _A  =  1  and  cp  =  0.7054.  Table  III  is  a  listing  of 


ML  at.  \  9\  o\  ■ 


-  b6  - 


P’a 

— 1  and  9'a  for  9  a  -  characteristics  at  A  where  ££pA#  =  5°  and  is 

tr  ** 


incremented  5° 

for  each  successive 

a  ~  characteristic 

0 

Table 

III 

i 

V, 

1 

P  ' 

A  i 

k 

i 

1 

* O8727 

I.I7I+5 

0*8727 

2 

.1745 

1.3491 

0*9599 

3 

"*2618 

1*5256 

I.OU72 

b 

.3491 

1*6981 

1.1345 

5 

.4363 

1.8727 

1*2217 

6 

.5236 

2.0^72 

1,5090 

7 

.6109 

2*2217 

1 • 3963 

8 

*6981 

2.3963 

1*1+855 

9 

.7854 

2,5708 

1*5708 

With  reference  to  Figure  6S  consider  known  point  Q(r^,z^) 
on  a  |B  -  characteristic  and  known  point  P(rp,zp)  on  an  a  -  characteristic 
Also  both  p  and  9  are  assumed  known  at  P  and  Q.  To  calculate 
R(r,z)j,  the  point  of  intersection  of  the  3  -  characteristic  through  P 
and  the  a  -  characteristic  through  Qs  the  following  procedure  was 
adopted*  A  first  approximation  of  the  actual  location  of  R(r9z)  is 
made  by  R^(r^sz^)  which  is  the  point  of  intersection  of  the  tangent 
to  the  3  «  line  through  P  and  the  tangent  to  the  a  -  line  through  Q* 

Both  r^  and  z^  are  calculable  from  the  simultaneous  equations 


;1  -  zp  =  (rp  -  rx)  cot  Cpp 
;1  -  ZQ  =  <rl  -  rQ}  tan 


and 


10 
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Intersection  of  Characteristics  Through 
Neighbouring  Points  to  Show  Approximation  for  True  Intersection  Re 


The  state  of  stress  at  R(rsz)  is  governed  by  (3o3<,8)  and 
Expressed  in  finite  difference  form  these  are 


k  (pi  •  P(p  +  2(tpi  -  cp0)  =  -2 


(rl‘r 0+Zl-Z q) 


ri+rQ 


and 


£  (Pp  -  Pp)  -  2(q>1  -  cpp)  =  -2  — t; 


(rl"rp°zl+zp) 


l+rP 


(4.4.5) 


(4.4.6) 


respectively,,  An  approximation  for  cp  at  R  is  made  by  cp^  at 
solved  from  (4o4<,5)  and  (404o6)„  A  closer  approximation  to  R  is  made 
by  niz^)  which  is  the  point  of  intersection  of  the  straight  line 

through  Q  whose  slope  is  the  mean  of  the  slopes  of  the  a  lines  at 
R^  and  Q  and  the  straight  line  through  P  whose  slope  is  the  mean  of 
the  slopes  of  the  (3  -  lines  through  R^  and  P„  Both  r0  and  z^  are 
determined  from  the  equations 


rfftUQldJ 
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cot  cp  +  cot  cp 

r2  -  rQ  =  (z2  -  ZQ)  — ' — “  (4,4.7) 

cot  cp  +  cot  cp 

and  z2  -  zp  =  (rp  -  r2) - - - ~  (4.4.8) 

A  second  approximation  to  cp  at  R  is  made  by  cp  solved  from  (4.4.5) 
and  (4.4.6) .  This  procedure  is  continued  until  the  difference  between 
two  successive  approximations  of  cp  is  less  than  the  nominal  accuracy 
of  the  calculations .  The  value  of  p  corresponding  to  the  final  value 
of  cp  is  then  derived  from  (4.4.5) » 

The  above  procedure  was  programmed  in  Fortran  for  the  I.B.M. 

1620  digital  computer.  The  p  and  q>  values  for  the  9  a  -  character¬ 
istics  at  the  singular  point  A  (Table  III)  and  for  the  14  points  along 
the  a  -  characteristic  AC  (Table  XX)  were  used  to  determine  the  fan¬ 
shaped  plastic  stress  field .  The  program  together  with  the  values  of 
p,  Cp,  r  and  z  at  each  point  computed  are  found  in  Appendix  I„  Figure 
3,  however,  contains  only  that  portion  of  the  field  that  is  applicable 
to  the  indentation  problem  considered. 

4.5  Plastic  Stress  Field  ADE . 

The  boundary  condition  of  zero  shearing  stress  along  the  straight- 
line,  smooth  boundary  AE  of  region  ADE  (Fig.  3  or  4)  causes  the  a 
and  (B  -  characteristics  to  meet  this  boundary  at  45°  angles.  Together  with 
the  a  -  characteristic  AD,  this  boundary  condition  determines  the  plastic 
stress  field  in  ADE .  This  arrangement  is  analogous  to  the  third  boundary- 
value  problem  [Hill,  1950  ]  in  plane  strain  plasticity. 

With  reference  to  Figure  4,  consider  point  T  on  a  -  line  AD 
at  which  p  and  Cp  are  known.  J(r,z),  which  is  the  point  of  intersection 
of  the  0  -  line  through  T  and  the  straight-line  boundary  is  sufficiently 


,  . 


I 
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located  in  the  following  manner 0  If  T  is  in  the  immediate  vicinity  of 
A,  then  an  adequate  approximation  to  J  is  made  by  J^(r^,z^),  the  point 
of  intersection  of  the  boundary  AE  and  the  straight  line  through  T 
whose  slope  is  the  mean  of  the  slopes  of  the  (3  -  lines  at  T  and  J„ 

Both  r^  and  zn  are  determined  from  the  equations 

Z1  =  R  -  r! 

(it. 5.1) 

and  =  cot 

where  R  is  the  radius  OA  of  the  cavity  at  surface.  The  pressure  p 
at  is  determined  by  (4,4o6)„  The  a  -  characteristic  J^K  is  then 

determined  from  the  data  on  a  ~  characteristic  AD  and  at  by  the 

procedure  described  in  4, 4*  The  remainder  of  the  field  in  ADE  is 
determined  by  repeating  the  above  procedure  upon  successive  replacement 
of  the  role  of  AD  by  J^K  „ 

The  Fortran  program  for  the  procedure  was  run  on  an  I.B.M.  1620 
digital  computer  using  the  data  for  the  14  points  on  AD  obtained  in  4,4, 
Originally  the  field  was  overdetermined  and  it  also  became  evident  from 
analysis  of  the  data  that  the  arc  length  corresponding  to  the  role  of  TA 
was  too  large  at  some  stage  in  the  computation.  This  rendered  the  approx¬ 
imation  of  the  point  on  AE  symbolized  by  J  by  the  point  symbolized 
by  inadequate,  A  more  dense  distribution  of  points  was  then  taken 

along  the  7th  a  -  line  along  AE,  To  accomplish  this,  large  scale 
graphs  were  constructed  of  the  variation  of  p,  cp  and  z  versus  r  along 
this  a  -  line.  From  these  graphs^  1 7  points  were  taken  and  another  program 
was  run  on  the  digital  computer.  Doing  this  served  the  dual  purpose  of 
retaining  the  accuracy  of  the  data  and^also^in  aiding  the  determination 
of  the  p  -  characteristic  EN  which  passed  through  apex  point  E  ,  This 


3  "  characteristic  EN  was  determined  by  interpolation,,  Graphs  of  the 
variation  of  cp  and  z  versus  r  along  MN  produced  were  constructed,, 
The  location  of  N  was  determined  by  trial  and  error  with  use  of  equation 
(Uc5cl)o  The  value  of  p  at  E  was  then  calculated  as  10.210  k  using 
(4„4o6)o  Figure  3  shows  only  that  portion  of  the  plastic  stress  field 
computed  from  the  data  that  is  applicable  to  the  indentation  problem,, 

The  Fortran  program  and  data  for  the  above  procedure  are  found  in 
Appendix  L 

b.6  Calculation  of  Indentation  Stress  from  Slip-Line  Field  of  the 
Incomplete  Solution. 

Figure  3  is  the  plastic  stress  field  or  equivalently  the 

slip-line  field  constructed  from  the  data  of  Appendix  I  and  represents 

the  incipient  plastic  flow  state  of  the  material.  Shown  also  in  Figure 

3  are  values  of  the  normal  stress  cr  in  units  of  k  at  13  points 

along  AE.  These  represent  the  distribution  of  the  normal  stresses  that 

must  be  supplied  by  the  conical  indenter  to  produce  this  incipient  plastic 

flow  state.  The  uniform  stress  <7  to  be  applied  to  the  conical  indenter 

z 

in  the  axial  direction  and  to  be  distributed  over  its  base  of  radius  R 
is  referred  to  in  this  thesis  as  the  indentation  stress. 


Reference  Diagram  for  Computation  of  the  Indentation  Stress 


, 


, 


With  reference  to  Figure  7>  the  indentation  stress  0"  is  given 


by 


2k 

R2 


f 


R 


J 


O'  r  dr 

n 


which  for  the  specific  value  of  R  =  25  becomes 


2  k 
625 


r25 

J  y(^)  ^ 


(4.6.1) 


where  y(r)  =  an  r  .  The  definite  integral  in  (4.6.1)  was  evaluated 
using  the  Gauss  quadratic  formula  (National  Physical  Laboratory  [ 1 96 1 ] ) . 
This  formula  states  that 


r 


j 


y(r)  dr  =  —  (b-a)  f  y(X)  dX 
^  \J 

-1 


n 


|  (b-a)  Y 


where  r  =  ~  (b+a)  +  —  (b~a)  X 


r=l 


(n) 


w^n)  y(x<">)  , 


are  zeros  of  Legendre  polynomials 


.(«) 


and  w  are  weights  such  that  the  formula  is  exact  when  y  is  any 

/  \  /  \ 

polynomial  of  degree  less  than  2n.  Values  for  X^.  '  and  w^  '  were 
available  for  n  =  16  from  the  University  of  Alberta  Computing  Center 
However,  the  points  along  AE  (Fig.  3)  at  which  or  is  known  were  not 
equally  spaced  nor  did  their  abscissa  values  coincide  under  the  trans¬ 
formation  X 


-  1  with  the  16  zeros  X^  '  of  the  Legendre  polynomials. 


This  thus  necessitated  interpolation  of  the  values  y(r)  at  points  having 
/  \ 

abscissas  X^.  under  this  transformation.  Interpolation  with  unequal 
intervals  required  the  calculation  of  some  Lagrangian  coefficients  (Comrie 
[1959])*  These  calculations  together  with  the  evaluation  of  the  definite 
integral  in  (4.6.1)  are  found  in  Appendix  II. 


Results  of  the  above  computations  gave  a  value  of  4.6424  k  for 
the  indentation  stress  er  .  Consequently,  the  yield  point  load  L  that 

Z  J 

is,  the  load  necessary  to  produce  incipient  plastic  flow  in  axially 
symmetric  identation  by  a  conical  punch  of  semi-angle  45° >  is  given  by 

L  =  4.6424  itkR2 

where  k  is  the  maximum  shearing  stress  of  the  indented  material  and  R 
is  the  radius  of  the  conical  cavity  at  the  surface. 

4.7  Discussion  and  Future  Requirements. 

The  value  obtained  for  the  yield  point  load  has  as  yet  not  been 
proven  to  be  the  actual  yield  load.  For  on  the  basis  of  limit  analysis 
theorems,  the  plastic  stress  field  is  unsuited  for  the  establishment  of  a 
bound  on  the  actual  yield  point  load.  The  field  has  not  been  extended 
into  the  rigid  region  and,  hence,  not  shown  to  be  statically  admissible; 
nor  has  any  attempt  been  made  to  show  that  an  associated  kinematically 
admissible  velocity  field,  compatible  with  the  stress  field,  exists  in 
the  deforming  region. 

Shields  [19551s  in  his  investigations  on  axially  symmetric  in¬ 
dentation  by  a  cylindrical  punch,  has  detailed  a  procedure  for  extending 
the  plastic  stress  field  into  the  rigid  zone.  This  was  done  using  a 
method  which  is  an  application  of  the  theorems  of  Bishop  [19531  on  the 
admissibility  of  incomplete  solutions.  It  is  proposed  that  a  direct 
application  of  Shield's  procedure  to  the  present  problem  in  conical  indent¬ 
ation  should  result  in  the  extension  of  the  plastic  stress  field  into  the 
rigid  zone.  Moreover,  the  resulting  stress  free  surface  should  lie  wholly 
within  the  boundaries  of  the  indented  material,  thus  meeting  one  of  the 
requirements  of  Theorem  1  of  Bishop's.  By  this  theorem^if  for  a  fully 


■  *  rt  -  J 


plastic  continuation  of  the  stress  solution  of  a  fully  plastic  region 
which  is  in  equilibrium  with  an  applied  external  load  and  which  contains 
a  compatible  deformation  mode,  the  resulting  stress  free  surface  lies 
wholly  within  or  upon  the  external  boundary  of  the  material,  then  a 
complete  solution  exists  and  the  actual  yield-point  load  has  been  found. 

To  meet  the  other  requirement  of  the  theorem,  a  velocity  field 
must  be  found  in  the  deforming  region  compatible  with  the  plastic  stress 
field,  A  method  by  which  this  may  be  done  is  now  considered.  The  region 
ABODE  (Fig,  3)  is  deforming,  whereas,  the  region  above  it  is  rigid.  Thus 
there  is  a  possible  velocity  discontinuity  across  the  3-line  BODE,  In 
terms  of  the  velocity  resolutes  U  and  W  defined  by  equations  (3o3°l8), 
it  is  possible  to  show  that  the  boundary  conditions  to  be  satisfied  on 
this  3-line  is  U  =  W  =  0,  If  now  the  boundary  condition  along  the  non= 
characteristic  line  AE  can  be  established  when  the  punch  is  penetrating 
in  an  axial  direction  at  a  rate  w^,  say,  the  entire  velocity  field 
can  be  determined  by  a  method  analogous  to  that  type  in  plane  strain 
(type  ( iii) ,  Hill  [1950]), 

In  accord  with  the  above  discussion,  the  validity  of  the  value 
4,6U24itkR2  for  the  yield  point  load  thus  rests  on  the  determination  of 
the  kinematically  admissible  velocity  field.  Until  this  has  been  done, 
this  value  is  only  a  tentative  one. 
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APPENDIX  I 


COMPUTER  PROGRAMS  AND  OUTPUT  DATA  FOR  PLASTIC  STRESS  FIELDS 


The  calculations  for  the  individual  plastic  stress  fields 
which  form  the  region  ABCDE  (Fig.  3)  were  greatly  facilitated  by  the 
use  of  an  I.B.M.  1620  digital  computer.  The  iterative  processes  by 
which  these  calculations  had  to  be  made  were  ideally  suited  as  computer 
problems .  The  Fortran  source  programs  for  each  field  were  designed  so 
as  to  be  adaptable  to  other  axially  symmetric  indentation  problems. 


A.  Plastic  Stress  Field  AC^D 
1.  Procedure 

As  detailed  in  4.4,  Chapter  IV,  the  actual  location  of  R(r,z) 
in  Figure  6  is  approximated  firstly  by  R^(r^,z^).  The  coordinates  of 
this  latter  point  are 


zi  - 


^rP*r<P  +  ZQ  COt  COt  ‘Pp  +  Z] 

cot  CpQ  cot  CPp  +  1 


and 


rl  = 


(Z]_  "  ZQ)  COt  +  rQ  » 


(A. 1.1) 


as  derived  from  equations  (4.4.4).  Now  the  state  of  stress  at  R(r,z) 
is  governed  by  the  finite  difference  equations  (4. 5)  and  (4.4.6). 
Solving  these  equations  for  cp  yields 


cp  = 


P„-P„  1  |lr-z)-(rp-zp)  (r+z)-(rn+zn) 


'p_12  +  I 

k  2 


r  +  r. 


r  +  r 


Lo  J  • 


cp  cannot  be  determined  explicitly  from  equation  (A. 1.2)  as  r  and  z 
are  as  yet  unknown.  However^  an  approximation  to  cp  is  attained  by 
say  cp^  upon  giving  the  values  of  r^  and  z^  to  r  and  z  respect¬ 
ively.  Then  a  second  approximation  to  R(r,z)  is  made  by 
whose  coordinates  are  given  by 


V' 

limb 


' 


CO  tCp^+CO  tCp  ^ 


r> 

c. ; 


1  z  +  (  cotep  +COtCp,  'Z 


\ 


„  4»  r 

1  /  p  p 


-  r 


and 


r_  = 


(Z2  "  ZQ 


cotffl^+cotq). 
_ >•;  . •>  + 


cottp^+cotcp^ 


r\ 

£L 


CO  tCPp4-CO  tCp  ^ 


(a.i.3) 


2  J  +  rQ  > 


as  derived  from  equations  ( 14- 0  4 <, 7 )  and  (4»4o8)  0  Upon  substitution  of 
these  values  of  r0  and  z^  for  r  and  z  in  (A. 1,2),  a  second 
approximation  to  cp  is  attained*  This  procedure  is  continued  until 
the  difference  between  two  successive  approximations  of  cp  is  less  than 
the  nomial  accuracy  of  the  calculations „  The  coordinates  of  R(r,z)  are 
then  taken  as  the  final  values  of  (A.l.j)  which  attained  this  condition. 
Finally ?  p  at  R(r,z)  is  solved  using  equation  (4<,4,5)° 


The  above  procedure  is  concerned  only  with  the  location  and 
state  of  stress  of  one  point  R(rj,z),  However,  it  is  possible  to  deter¬ 
mine  an  entire  (B-line  (or  a-line)  easily  with  a  computer  by  transferring 
the  role  of  points  and  introducing  others.  For  at  singular  point  A 
(Fig,  4)?which  is  to  be  considered  as  a  degenerate  (3- line5  there  is  init¬ 
ially  the  first  of  9  points  (Table  III)  to  play  the  role  of  Q  and  each 
of  14  points  (Table  IX)  along  AC^  to  play  the  role  of  P,  After  R(r,z) 
has  been  calculated^ it  is  transferred  to  act  as  P  and  the  next  point 
at  A  is  introduced  to  act  as  Q,  In  this  manner  the  first  a-line  in 
AC pD  is  determined.  The  arrangement  then  becomes  similar  to  the  orig¬ 
inal  one  with  the  newly  computed  a-line  acting  as  AC^,  By  repeating 
the  entire  procedure  the  stress  field  is  determined. 


2,  Fortran  Source  Program,  Input  and  Output  Data, 

The  equations  of  A, 1  were  translated  into  Fortran  language  of 
the  computer  using  the  following  notation  for  the  various  quantities 


,  - 


% 


appearing  in  them” 


X(I)  -  rQ  ,  Y(I)  =  20  , 

T(I)  =  <P0  , 

P(I)  =  p 

XX  =  rp  ,  Ti  =  zp  , 

XXX  =  rp  ,  YYY  = 

XXXX  =  r0  ,  YYYY  =  z2 

TT  =  cpp  , 

pp  =  Pp 

w  =  cot  CpQ  ,  V  =  cot  cpp 

b  =  , 

c  =  (r+z)-(rp 

+zp) 

r+rQ 

r+rp 

TA  =  COtVC°tCpl  TR  2 

1  .D  =  "  “  — -  9 

2  eotcpp+cotcp^ 


The  Fortran  source  program  ^which  embodies  the  procedure  of 

A. l^is  found  in  the  following  attached  sheet  together  with  the  input 
data  which  gives  the  locations  and  state  of  stress  at  the  known  points 
along  AC^  and  the  singular  point  A.  These  are  followed  by  the  output 
data  which  determines  the  fan-shaped  plastic  stress  field  AC-^D. 

B.  Plastic  Stress  Field  ADE. 

For  the  remaining  region  ADE  (Fig.  4)  the  procedure  adopted 
was  essentially  the  same  as  the  one  used  above  except  that  successive 
a=lines  rather  than  0-lines  were  computed  to  extend  the  stress  field 
AC^D  into  the  region.  For  this  reason  detailed  discussion  of  the  pro= 
cedure  would  be  unwarranted.  Thus  only  the  Fortran  source  statements 
and  the  input  and  output  data  for  ADE  are  to  be  found  in  the  following 


attached  sheets. 


- 

1 

Fortran  Source  Statement  for  Plastic  Stress  Field  AC^D 


n  i  mfns  t  on  x  (  ^ »  y  (  ^  n  ) ,  t  (  ^  ,  p  ( c  p  ) 

2  RFAH*N*M»JJ*F 

nO  i  I  =  l  *  N 

1  READ*X(  I)  *Y(  I  )  *T(  T)  *P(  T  ) 

00  ?0  J  =  J  J  » M 
REAn,XX»YY»TT *PP 
DO  A  0  1=1  *N 
V  =  G  .  3 

W=COS (T(I))/SIN(T(I)) 

U  =  C  0  S ( TT) /SIN(TT) 

YYY= ( ( ( XX~X( I ) )+Y ( I )*W) *U  +  YY ) / (U*W+1 .0 ) 
XX X= ( YYY-Y ( I ) ) *W  +  X ( I ) 

5  R= ( XXX~VYY-XX+PY ) / ( X XV+XX ) 

C  = ( XXX+YYY-X (  I  )  — Y (  I  )  )  / ( XXX+X (  I ) ) 

P  H  I  = ( T (  I  )+TT)/?.0-(PP-P(  I  )  ) / 4 • 0  + ( B  ■  ")/?. 
F=  ARS ( V-PH I  ) 

IF  (F-F)7,7,R 
R  V=  PH  I 

T  A= ( W+COS ( V ) /S I N ( V )  ) / 2 • 0 

TB  =  ?  •  0/ ( U  +  COS ( V ) / S I N ( V )  ) 

YYYY=(TA*Y( I ) +TB*YY+XX“X ( I ) ) / ( TA+TB ) 

XX XX= ( Y YYY-Y ( I ) ) *TA+X ( I ) 

yyy=yyyy 

xxx=xxxx 

GO  TO  R 

-7  pf=  (  T  (  T  ) -r-PHI  )  *?  .  3+P  (  T  ) 

IF  ( S c N S F  SWITPH  8)9, 1  0 

9  PRINT , I ,J*PF,PHI ,XXX,YYY 

1 0  PUNCH,  I  ,  J  *PF ,PHT 
PUNCH ,  I  ,  J  ,  XXX , YYY 
X(  I  )=XXX 

Y  (  I  )  =  Y  Y  Y 

T (  I  )=PHI 
P (  I  )  =PF 
XX=XXX 
YY  =  YYY 
T  T  =  PH  I 
oQ  PP  =PF 

20  rONUNUF 

PAUSF 
GO  TO  2 
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r  .  O  (  r  H 

- - - —  „T..  r-;  .  .  ’  ,  j  '  7-;.t  j  q  q - p- 

THQ. 

- f  XT'  •  It  »  t^Ot*  'H - 

— Y-^=-f  I  )  Y - 

- T77  =  YY - 

T  H^=  T  T 
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Stress  Field  AC-^D 
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Input  Data  for  Plastic 


9  14  1  0.00001 


Data  at  Singular  Point  A 


2  5. 0000 

0  .  "i  o  00 

0 . p  727 

1 

.  1  745 

25.0000 

0 . 0  0 

C  .  Q5QQ 

] 

.0491 

2  5.0  0 

0  .  o  n  o  0 

i , 0472 

1 

.  *  286 

25.0000 

0.0000 

1  .  1  ?  4  5 

1 

.  OOP  1 

25.0000 

0.0000 

1.2217 

1 

.  8727 

25. 

0. 

1 .5090 

2 

.  9472 

2  5.  00 

0.0000 

1.3963 

2 

.2217 

2  5. 

0.1 

1.4835 

2 

.9963 

25.0000 

0  . 

1 .5708 

2 

.  8708 

Data  for  a-character istic  AC^ 


25.5190 

0 . 5 1  0  5 

C  .7654 

0.9992 

26.0866 

1  .  04  40 

0.74^9 

0.9968 

26.7068 

1  .  6042 

0.77^0 

0.992  7 

27. 3896 

2.  1957 

0.7044 

0. 0868 

28.1436 

2.8237 

0 .6833 

C  .  9  7  0  1 

18.9876 

3.4952 

0.66  17 

0.9693 

29.9393 

4.2190 

0.6393 

.  '0  7  2 

3  1.0262 

5.0070 

•  6  1  6  0 

0.9423 

32.2847 

5.8749 

0 .5914 

0.9243 

33.7732 

6.8462 

0.5652 

0.9023 

35.5813 

7.9567 

0.5367 

.  '  5  i 

37.8677 

9.2667 

r .5048 

0.8408 

40.9459 

10.8925 
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OUTPUT  DATA 


(3 -line 


i 

p=  l*i  /  3  5  0  j  1 

(p=  .851818 

1 

r=s  25*4  4 

Z—  .  3  g  0  J  2  1 

1 

1.351  .1 

. 

1 

.42 

.5 

1 

1.3283398 

1 • 0248496 

1 

25. 

.62886672 

1 

1  •  7  0  3  26  94 

1 .  x  1x5096 

1 

25.316862 

•  6  3  0  5  6  3  4  3 

1 

1*8620332 

1.1982410 

1 

Z5 . 258646 

.68533327 

1 

^;*\j^89  /  7  8 

1.2850972 

1 

2  3  *  1982  7  z 

•  703  9  9  7  90 

1 

2.2359746 

1 . j  7205  77 

1 

25.136z23 

•72135883 

1 

2.4129740 

L  0  0  c 

1 

25.073010 

.7312691 

1 

2 . 5  9  c  u  6  0  6 

1.54627 5 

1 

25.009018 

.  564251 

2 

1.1747827 

•  8  3  0  7  72  45 

2 

25.994892 

1.13! 5563 

2 

1.3530607 

•916 36 58^ 

2 

23 . 89491z 

1.2195363 

2 

1  •  551682b 

1  •  002  14^+4 

2 

2  5.  (87724 

1 .294925c 

2 

1.  7 1 0  5  7 1 2 

1.0880969 

2 

25.674105 

1.3609993 

2 

1.8896474 

1.1742110 

2 

25.554897 

1.4171  ;h7 

2 

2.0689856 

i  .260j4  7  0 

2 

25.430819 

1.4630646 

2 

2.248 5^24 

1 .34  709  c6 

2 

25.302 8 0 

1.4981463 

2 

2.4281108 

1 .4358264 

2 

25.171857 

1  •  5  2  2  0  6  7  5 

2 

2.60/8/94 

x . vlidl  I  7 

2 

25.038 <  9  / 

1.5 34365  7 

3 

1.1722828 

. Ocl454lZ 

3 

26 .3674  78 

1 • 7492149 

3 

1 .3523402 

• 

3 

zo .41 5558 

1 . d  bxb i  71 

3 

1 .5329180 

•  9  7  8  993  06 

3 

26. 252088 

2.0032075 

3 

1 .7139326 

1  •  w  6  41 3  74 

3 

26.07816^ 

2.10  9248 

3 

1  .8952 >6z 

1 . 1 4  5  5  2  6  5 

3 

2  5 . 8  9  4  7  2 

z.zOib^w j 

3 

2.0/70/6^ 

1 • z  3  5  22  7  3 

3 

25.70357 9 

<i.<x  7  80  5  5  9 

3 

2.2591758 

1 . 5  Z  1  Z  4  5  0 

3 

23.505346 

<- .  3  3 

3 

2.4414956 

1.40753x4 

3 

23 . 3u 1 72  3 

x.  .  38060x4 

3 

2.6240904 

1.4^4xz > 1 

3 

23 • 09400O 

4 

1 .1678951 

. 78775593 

4 

27.2^1459 

2.4jc5190 

4 

1.3496593 

•8 7 136116 

4 

2o. ^96080 

X.500U1JO 

4 
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FORTRAN  SOURCE  STATEMENT  FOR  PLASTIC  STRESS  FIELD  ADE 


D I MFNS ION  X ( 50 ) * Y ( 50 )  ,  P { 50 )  ,  T ( 50 ) 

2  REAO,N»M»JJ,K*F 
DO  1  1  =  1  ,  M 

REAO,P (  I  )  ,  T (  I  ) 

1  RFAO,X ( I ) , Y ( I ) 

00  20  J  =  J  J  * N 
T  T  =  T  (  1  ) 

T  2  =  1  .  5708 

T (  1  )  =  (T  (1  )  +T  2 ) / 2 • 0 

17  U  =  fOS( T ( 1  )  ) / S I N ( T ( 1 )  ) 

X  2  =  (  Y  (  1  )  +  X  (  1  )  *U-2  5. 0  )  /  (  U-l  .‘0  ) 

Y2=?5.0-X2 
]8  T ( 1 ) = T T 

R= ( X2-V2-X ( ]  )  +  8 ( 1 )  ) *2.0/ ( X2  +  X ( 1 )  ) 

P2= ( T 2-T ( i ) ) *2 . 0  +  P ( 1 ) -7 
IF  (SENSE  SWITCH  3)6,7 

6  PR INT* J,K,X2*Y2>P2*T2 

7  PUNCH,  J,K,P2»T2 
PUNCH, J ,K , X2 , Y  2 
no  15  1=2, M 
V=0. 0 

S  =  COS ( T2 ) /SIN( T2 ) 

R  =  COS ( T ( I  )  )  /  SIN ( T (  T  )  ) 

Y  8  = (  ( X (  T)-X2  +  Y2*S)*R  +  Y(  I  )  )/( S*R+i . 0 ) 

X3= ( Y3-Y2 ) *S+X2 

10  D= ( X8-Y8-X (  I  ) +8 (  I  )  ) / ( X^  +  X (  I  )  ) 

C= (X8+Y8-X2-Y?)/(X7+X2) 

T  8  = ( T 2+  T (  I  )  ) / 2 • 0— ( P (  I  )-P2 ) /4.0+ ( O-C ) /2 .0 
F= APS ( V-T  3 ) 

I F  ( F-F  )  8 ,8,9 
9  V  =  T  3 

R A= ( S  +  COS ( V) /SIN( V) ) /? . 0 
R  r  =  2 . 0 / ( R  +  COS ( V ) / S I N ( V )  ) 

Y  4  =  (  Ra*Y2  +  Rp*Y  (  I  )  +  X  (  I  )-X2)/(RA  +  RD) 

X4  = ( Y4-Y? ) *  P  A  +  X  2 

Y3=Y4 
X  8  =  X  4 
GO  TO  10 

r  P3= ( T2-C-T3 ) *2. 0+P2 

IF  (SENSE  SWITCH  3)11,12 

11  P  R  I  N  T  »  J , I  »  X  3 ,  *  "  ,  P  3 ,  i  -  , 

12  PUNCH, J , I , X3 , Y3 
PUNCH, U, I  , P  3 , T  3 
X2  =  X3 

Y2=Y3 
P2  =P7 
T  2  =  T  3 
X (  T - 1 ) =X3 
Y (  1-1  ) =  Y  3 
P  (  1-1  ) =  P  3 
1  5  J  (  j  - 1  )  =  T  3 
20  CONTI NUF 

2 

END 
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INPUT  DATA  FOR  PLASTIC  STRESS  FIELD  ADE 
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2.6614754 

10 
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OUTPUT  DATA  FOR  PLASTIC  STRESS  FIELD  ADE 
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26,893636 
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1 .0118 2 8 9 

1 1 
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2.9169328 
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17.385787 

3.4068618 
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. 80725784 

10 
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FORTRAN  SOURCE  STATEMENT  FOR  PLASTIC  STRESS  FIELD  ADE 

(Run  2) 

0  I MFNSION  X(50) >Y (50) , P ( 50  ) »  T (50) 

2  READ*N*M* JJ»K*F 
DO  1  1  =  1  » M 

RF AD, P (  I  )  , T (  I  ) 

1  RFAD,X(  I  )  ,Y(  I  ) 

DO  ?_0  J  =  J  J  *  N 
T  T  =  T ( 1  ) 

T2  =  i  . c  7 CP 

T (  1  )  =  ( T  (  i  )  +  T2 ) /? • 0 
1 7  U  =  COS( T (1  )  ) /SlN( T ( 1  )  ) 

X  2  =  ( Y ( 1  )  +  X ( 1  ) *U~2  5 • 0 ) / ( U-l  . 0 ) 

Y  2  =  2  5  •  0  -  X  2 
1 p  T (  1  )  =TT 

B= ( X2-Y2-X ( 1  )+Y(]  )  ) *  2 . 0 / ( X  2  +  X ( 1  )  ) 

P  2  = ( T  2~  T ( 1  )  ) *2 . 0  +  P ( 1 ) -B 
IF  ( SFNSF  SWITCH  3)6*7 

6  PRINT  »J»K>X2»Y2>P2>T2 

7  PUNCH , J ,K » P2 * T2 
PUNCH »  J  *  K  »  X2  »  Y2 
DO  l c  I  =2  * M 
V=0.0 

S=C OS ( T  2 ) / S I N ( T  2 ) 

R  = COS ( T (  I  )  )  /S  I  N ( T (  I  )  ) 

Y  p  = (  ( X (  I ) -X2  +  Y2*S ) *R  +  Y(  I  )  ) / ( S*R+i . 0 ) 

X  2  = ( Y2-Y2 ) *S  +  X2 

10  D= ( X3~Y7~X (  I  ) +Y (  I )  ) / ( X2  +  X (  I  )  ) 

C  =  (X2  +  Y7-X2-Y2)/(X'5  +  X2) 

T  7  =  ( T  2  +  T (  I  )  )  / 2 . 0— ( P(  I  )-P2) /4.0+( D-C ) /2 .0 
F= ABS ( V-T 2 ) 

I F  ( F-F  )  0  »  «  ,9 
O  V=  T  2 

R  A  = ( S  +  rOS ( V) /SIN ( V)  ) / 2 . 0 
RB=2.0/(R+COS(V)/SIN(V) ) 

Y4= ( RA*Y2+RB*Y ( I )+X ( I ) -X2 ) / ( RA+RB ) 

X4  = ( Y4-Y2 ) *R A  +  X2 

Y  2  =  Y4 
X3=X4 

GO  TO  10 

p  P3= ( T2-C-T3 ) *2 . 0+P2 

I F  ( SFNSE  SW I TCH  3)11,12 

11  PR  INT,J,I  ,X2,Y2,P3,T3, 

1 2  PUNTH, J, I >X7 ,  Y7 
PUNCH  » J  *  I  ,P3»T2 
X  2  =  X  2 

Y  2  =  Y  2 
P2  =  P2 
T  2  =  T  2 

X  (  I-]  ) =  X2 

Y  (  I-l  ) =  Y3 
P  (  1-1  ) =P2 

13  T (  I  —  1  ) =  T  3 
20  CONTINUE 

P  A  U  S  F 
GO  TO  2 
c  N  D 
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INPUT  DATA  FOR  PLASTIC  STRESS  FIELD  ADE  (Run  2) 


1  7  17 

Q. 0  0  0  0 1 


p=3 .5044 
r=l  7  •  2  5  0  0 
3 .5022 

1 7.2775 

3 .4964 
1  7  .  ^00 
3 . 4372 
17.4626 
3.4724 
17.650 
3 . 4572 
17.8431 
3.4453 
18.0000 
3 . 4301 
1 8.2000 
3.4154 

18.4000 
3 . 4069 
18.5168 
3 . 3900 
18.7500 
3.3723 
19. 0  u j 0 
3.358 
19.2. 

3 . 3450 

19.400 0 
3.3385 

19.5 

3.33q4 

19.6000 

3 . 3284 

19.6  r  3 2 
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.  536  6  =cp 
0 • 02  00=z 
.5167 
9.6713- 

.4787 

1  0 . 6  4 

.4558 
1 1  .8^00 
.41  90 
13.2000 
.  3882 
14.3850 
.  3684 
] 5. 1 800 
.  34^0 
16.1 200 
.  3237 

16.9850 
.7127 
17.7858 
.  2  9  0  0 
18.2200 
.  2706 
1  °  ~  6 
,2560 
19.7200 

.2420 

2  0 . 3  8 
.7787 
20.6400 

.  7  2  8  ^ 
20.9250 
.2267 
2  1. 
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OUTPUT 

DATA  FOR  PLASTIC 

STRESS  FIELD 

ADE  (Run  2) 

a- line 

(3  -line 

r 

z 

P 

CP 

1 

1 

24.255116 

. 74488400 

2.6700839 

1  «  5  7 

1 

2 

24.26581 8 

1.5640013 

2. 6881425 

1 .5446685 

1 

3 

24.301383 

2.4572810 

2. 7045221 

1  .  5  1  7  3  c  8  P 

1 

4 

24.367902 

3.4366166 

2.7}p99q5 

1 .4886  3  07 

1 

5 

24.473202 

4.5143963 

2.7313561 

1,4582249 

1 

6 

2 4.628340 

5.7121557 

2. 7410509 

1 .42 c 8 245 

1 

7 

24.848038 

7.0527746 

2 • 74762  71 

1 ,3909999 

1 

8 

25.1 53235 

8 . 568r'  163 

2. 7500981 

1 . 3533466 

1 

9 

25.574217 

10.302651 

2. 7475369 

1*3121430 

1 

10 

26.157828 

12.320157 

2.7381103 

i  . •' 

1 

1 1 

26.979745 

14.719267 

2.7 1Q3391 

1.21 53447 

1 

12 

28.1  75567 

17.667646 

2.6871059 

1*1563 

] 

13 

30.018171 

21.480236 

2.6332539 

1,086071 7 

1 

14 

33.24912 2 

26.978370 

2. 5387225 

. 99636ro0 

2 

1 

23.425014 

1.5749860 

2.7761267 

1  .  <=  -7  n  o  o  o  0 

2 

2 

23.43798P 

2.4925362 

2. 7927407 

1,5426379 

2 

3 

23.481232 

3 . 4987907 

2.8073335 

1 . -129722 

2 

4 

23.562882 

4.6065980 

2.8196787 

1.481515 P 

2 

5 

23.694030 

5.8380428 

2. 829162c 

1  .44_7°40? 

2 

6 

23.889684 

7.2167072 

2. p  9^2771 

1 .4117077 

2 

7 

24.171155 

8.7755923 

2. 836955  5 

1 .3726663 

2 

8 

24.56934° 

10.559069 

2.8391623 

1,3298 

2 

9 

25.131978 

1 2 .633384 

2.8219147 

1  .  2  p  2  8  ^  0  0 

2 

10 

25.9361 34 

1 5.098581 

2 . 8005231 

1  .2290464 
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1  1 

27.119742 

18.125041 

2. 7645511 

1,1676809 
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12 

28.959766 

22.031740 

2.7053339 

1 #0948 1 6 

13 

32.206993 

2  7  •  6  5  C 

2.6027195 

1  .001 1 

14 

32  •  207372 

27.650690 

2.6026975 

1  •  0  0 1 1 

1 

22.494176 

2.5058240 

2.8907371 

1  , 

2 

2  2 . 5  1  0  1  2  1 

3.5418467 

2.9055395 

1  .  5400240 

3 

22.563954 

4.6830567 

2.9179181 

1.5073219 

4 

22.666964 

5.9521767 

2.9272015 

1 .4723 4 4 2 

5 

22.834547 

7.7735947 

2.9328039 

1 .434621 1 

6 

23.088530 

8.9811441 

2.9335563 

1 . 3937092 

7 

23.460652 

10.820416 

2.9282845 

i #348? 

8 

23.999803 

12.959016 

2.9148189 

1 . 2  0 ' 

Q 

24. 7849c0 

15.498869 

2. 8902071 

1  .2488000 

10 

25.957050 

18.613072 

2.8496265 

1  ,179117 9 

.1 

27.798299 

22.624391 

2. 7838143 

1 ,103 

12 

31  .071035 

28.373621 

2.6713677 

1 .0061 i 

13 

31.071434 

28.3742 c2 

2.671 3449 

1.0061141 

14 

31.072801 

28.37641 1 

2.671 2621 

1 .0060°P8 

1 

21.441714 

3.5582860 

3.0164567 

• 

2 

21  .461  877 

4.7369361 

3. 02QQ799 

i #536596 4 

3 

21.530flP3 

6.0486 1^0 

3.0381149 

1  .  4  0  0  g  ]  4  ^ 

4 

21 .664657 

7.5186734 

3.0431737 

1.460255* 

5 

21.885755 

9.1819078 

3.0428509 

1 .4171492 

6 

22.227056 

11.085290 

3. 03c7Q1 5 

1  . 
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22.739076 

1  3. 2978°7 

8.01^58  O'3 

i  .°i 
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2  3 . 5  0  3  3  1  3 

1*5.923  5  96 

2.9909297 

1 • 2  n  0  3241 

9 

24.664736 
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1  0 

26.51 2313 
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1 1 

29.822920 

2  9 . 1  6  0  9  9  c 

2. 7447621 

12 

29.823320 

29.161633 

2. 74c  738 1 

13 

29.824691 

29.163816 

2. 7456515 

14 

29.827P47 

29. 1 68P40 

2. 74^441 9 

1 

20.242206 

4.7577940 

3. 1 448790 

2 

20.268601 

6.11 79272 

3. 1 660130 

3 

20.360060 

7.6439472 

3. 1 704580 

4 

20.5401  p  9 

9.3718435 

3.1690206 

5 

20.8434^1 

1  1 .3501  89 

3  .  1  5  9  7  9  2 

6 

21.322725 

13.649677 

3. 1401304 

7 

22.062957 

16.376142 

3.1062970 

8 

23.214565 

19.707721 

7.0423794 

Q 

25.075403 

23.97296 2 

2*96754* ' 

10 

28.440739 

30.027203 

2. 82^3062 

1  1 

28.441140 

30.027846 

2.82 72806 

12 

28.442521 

30.030055 

2.8271890 

13 

28.444404 

70.034177 

2.  P  2 4 9 6 7 8 

14 

28.451976 

30.045178 

2.8265074 

1 

18.854647 

6.1453530 

3. 3172832 

2 

1 8.890642 

7.7767351 

7.7214834 

3 

19.0171 PP 

9.5412406 

7. 3] 76388 

4 

19.271139 

1 1 .609406 

3. 3067724 

5 

19.708741 

14.013726 

7. 2826504 

6 

20.42011 4 

16.861929 

7.2419144 

7 

21  .467274 

20.777926 

4. 1 77ft4l P 
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23.448045 

24.758962 

3.0785580 
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26.893233 
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2.01 77430 

1  .  0098890 

1  ,6698843 

1 .OOQR^^P 
1 . S7OOO0O 
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1 .4899232 

1  .  444041 3 
1.393  21 0 

1.3374532 
1 .2744658 
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1 , 1 1 7Q£ 
1,0120716 
1  .0120661 
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1  .  0 1 2  0 1 44 

1.0119577 

1 . 57OP000 
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1  .4060440 
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1  .34050^4 
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1  .21  427P1 
1 .1237374 

1  .  01  1  «344 
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10 
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APPENDIX  II 


.THE  EVALUATION  OF 


r  25 

I  a  r  dr  BY  NUMERICAL  METHODS 

J  n 

0 


For  the  determination  of  the  indentation  stress  cr  in  4.6, 


Chapter  IV,  the  value  of  the  definite  integral 


25 


a  r  dr  was 
n 


0 


required.  This  integral  was  evaluated  using  the  Gauss  quadrature 
formula  (National  Physical  Laboratory  [1961]). 


By  this  formula 

r 

0 


a  r  dr 
n 


J 

0 


r  25 

I  y(r)  dr 


-  s§  f\ w 

^  J 

-1 


dX 


n 


=  7  y(x(n))  , 


r=l 

where  y(r)  =  cr^  r  ,  X  =  “  -  1  ,  X^  '  are  zeros  of  Legendre  poly¬ 
nomials  and  w£n^  are  weights.  Values  of  X^n^  and  w^n^  corresponding 
to  n  =  16  that  were  used  in  this  formula  are  listed  as  follows: 

(n) 


X' 


n 


w 


“0. 989^0 

0.02715 

“0. 94458 

0.06225 

-0.86563 

0.09516 

-0.755^0 

0.12463 

-0.61788 

0.14960 

-0.45802 

O.I69I6 

-0.28160 

0.18260 

(continued) 


-  81  - 


x(n> 

n 

w 

r 

r 

-0.09501 

0.18945 

0.09501 

0.18945 

0.28160 

0. 18260 

0.45802 

0.16916 

O.6I788 

0.14960 

0.75540 

0.12463 

O.86563 

0.09516 

0.94458 

0.06225 

0.98940 

O.02715 

AlsOjthe  values  of  r,  cr^,  X  and  y(X)  for  22  points  along  AE  (Fig.  3) 
are  given  as  follows  where  the  value  of  unity  has  been  assigned  to  k? 


r 


<r 

n 


x  y(x) 


0 

11.210 

-1.00000 

0 

0.5926 

8.8907 

-0.95255 

5.2955 

2.2422 

7 • 4639 

-0.82052 

16.7556 

3.9224 

6.7594 

-0.68621 

26.5151 

5-5884 

6.5549 

-0.56895 

34.1350 

6.0957 

6. 1657 

-0.51254 

57.5721 

7 • 5707 

5.8955 

-0.41054 

45.4540 

8.7000 

5.6552 

-0.50400 

49.2002 

9.7621 

5.4851 

-0.21905 

55.5266 

11.2590 

5.2608 

-0.09928 

59.2515 

12.8661 

5.0405 

0.02929 

64.8516 

14.2195 

4.8742 

0.15754 

69.5077 

15.2719 

4.7580 

0.22175 

72.5585 

15.9578 

4.6524 

0.27662 

74.2421 

(continued) 


r 

er 

n 

X 

y(x) 

17.2256 

4.5062 

0.57805 

77.6220 

18.8546 

4.5175 

0.50857 

81.4010 

20.2422 

4.1569 

0.61958 

84.1448 

21.4417 

4. 0165 

0.7155^ 

86. 1206 

22.4942 

5.8907 

0.79954 

87.5182 

25.4250 

5.7761 

0.87400 

88.4551 

24.2551 

5.6701 

0.94041 

89. 0186 

25.0000 

5.5708 

1.00000 

89.2700 

Since  no  values  of  X  listed  above  coincided  with  any  value  of  X^n^, 

(n) 

it  was  necessary  to  interpolate  the  values  of  y(X^.  ').  This  was 
accomplished  using  the  Lagrange's  method  of  interpolation  with  unequal 
intervals  (National  Physical  Laboratory  [I96I])  using  4  consecutive 
values  of  y(X)  in  any  one  interpolation  calculation.  The  Lagrange's 


method  thus  used  may  be  written  as 


y(X^n))  =  L(X_1)y(X_1)+L(X1)y(X1)+L(X2)y(X2)+L(X3)y(X5) 

or  any  variation  thereof  depending  upon  the  value  of  X^n^  used  in 
relation  to  the  4  values  of  X.  The  L's  are  the  Lagrange  coefficients 
and  extensive  tables  for  them  are  available  only  for  equal  interval 
interpolation.  However since  the  values  of  X  form  unequal  intervals, 
this  necessitated  the  computation  of  new  Lagrange  coefficients  by  a 
simple  technique  outlined  by  Comrie  [19591*  In  way  of  an  illustration 5 
the  calculation  of  y(X^)  is  given  here;  the  values  and  not  the  cal- 
culations  of  the  other  y(X/.n')  being  given  later. 


c^8.8l 

“  33  “ 


x 


i 

r 


0.9840 


X 

y(x) 

DIFF. 

NUM.xlO6 

c 

DEN.xlO 

*  L(X) 

L(X)y(X) 

X  =  -1.00000 

0 

0.01060 

1,897.06 

2,683.60 

+0.706909 

0 

x2  =  -0.95235 

5,2953 

0.03705 

542.75 

1,671.81 

+0.324648 

+1.7191 

x,  =  -0.82052 
3 

16.7356 

0.16888 

119.07 

3,177.89 

-0.037468 

-0.6270 

x^  =  -0.68621 

26.5131 

O.30319 

66.32 

11,216.51 

+0.005913 

+0.1568 

y(X^.)  =  +1.24089 


The  values  of  the  columns  not  se If “explanatory  were  evaluated  under  the 
following  scheme  as  given  by  Comrie. 


X  =  n 
r 


X 

DIFF. 

NUM. 

DEN. 

L(X) 

X1  - a 

n-a 

(b-n) (c-n) (d-n) 

(b-a) (c-a) (d-a) 

+ 

x2  =  b 

b-n 

(n-a)(c-n)(d-n) 

(a-b)(c~b)(d-b) 

+ 

Xj  «  c 

c-n 

(n-a) (b-n) (d-n) 

( a-c ) (b-c ) (d-c ) 

- 

X.  =  d 

4 

d-n 

(n-a) (b-n) ( c-n) 

(a-d)(b-d)(c-d) 

+ 

In  this 

example 

n  lies  between  a 

and  b,  but  the  method  is 

NUM . 
DSN. 


for  other  positions  of  n.  The  sign  of  L(X)  is  determined  by  assigning 
+  to  the  two  adjacent  values  of  n  and  then  alternating  the  sign  in 
each  direction. 


The  remaining  values  of  y(X^n  )  are  included  in  the  following 
tabulation  from  which  the  definite  integral  under  consideration  was 


evaluated . 


A 

J 


Therefore 


-  81+  - 


y(X^n)) 

(n) 

w 

r 

W(n)  y(X(n 
r  r 

1.2489 

0.02715 

0.03391 

6.0872 

0.06225 

0.37895 

13.0659 

0.09516 

1.24334 

21.6341 

0. 12463 

2.69621 

31.0420 

0.14960 

4.64376 

40.7511 

O.I69I6 

6.89333 

50.3632 

0.18260 

9. 19647 

59.4181 

0.18945 

11.25682 

67.6350 

0.18945 

12.81352 

74.4173 

0.18260 

13.58882 

80.0079 

0.16916 

13.53390 

84.1110 

0. 14960 

12.58267 

86.8268 

0. 12463 

10.82105 

88.3745 

0.09516 

8.40963 

89.0452 

0.06225 

5.54342 

89.2444 

0.02715 

2.42316 

£  W^n)  y(X^n))  =  116.05896 


2^  ~ _ , 

f  y(r)  dt  =  )  w^n  y(x^n')  =  Ht50.Y3T00 


4.61+24  k  . 
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